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ABSTRACT 


The  heating  of  ions  in  a  magnetized  plasma  by  the  use  of  the 
second  order  fields  generated  by  the  nonlinear  mixing  of  two  whistler 
modes  is  examined.  The  resulting  kinetic  equations  describing  the 
mixing  and  heating  process  are  solved  using  the  method  of  orbit 
integrations.  Two  techniques  are  available  for  the  optimization 
of  the  energy  dissipated  by  the  ions.  One  is  to  allow  the  mixed 
wave  to  approach  a  natural  mode  in  the  plasma,  resulting  in  a  field 
resonance.  A  second  method  is  to  have  the  ions  absorb  energy  through 
cyclotron  damping.  The  resulting  sensitivities  to  fluctuations  in 
frequency,  density,  static  magnetic  field  and  inclination  of  the 
incident  sources  are  given  for  the  two  techniques. 
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CHAPTER  1  INTRODUCTION 


The  present  state  of  research  indicates  that  the  most  promising 
heating  techniques  for  the  attainment  of  controlled  thermonuclear 
conditions  are  those  for  which  the  energy  is  selectively  coupled  to 
the  ions.  Techniques  dependent  upon  ion-electron  collisional  effects 
become  inefficient  for  plasmas  with  temperatures  greater  than  lOOev. 
For  this  reason,  heating  schemes  which  employ  an  ion  cyclotron 
resonance  effect  are  of  great  interest. 

In  this  chapter,  a  summary  of  the  techniques  that  have  been  used 
to  date  will  be  given,  followed  by  a  discussion  on  the  method  of 
difference  frequency  harmonic  ion  heating. 


1.1  Heating  Techniques  Used  to  Date 


Some  of  the  schemes  that  have  been  employed  to  date  in  the  heating 

of  plasmas  are  enumerated  below.  Those  involving  magnetic  pumping  have 

been  described  by  Berger"*'.  A  survey  of  experiments  in  which  an  ion 

cyclotron  resonance  effect  is  employed  for  ion  heating  is  given  by 

2 

Hooke  and  Rothman  . 

(i)  Ohmic  Heating.  The  use  of  ohmic  heating  in  the  stellarator 

3 

experiments  is  described  by  Rose  and  Clark  (pp.  440  to  449).  Plasma 
heating  arises  strictly  from  the  resistivity  of  a  plasma  to  a  current 
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which  is  produced  by  a  unidirectional  electric  field  pulse.  This 

method  of  heating  will  be  efficient  only  at  low  temperatures,  because 

-3/2 

the  plasma  resistivity  is  proportional  to  T  ,  where  T  is  the 
temperature  of  the  plasma.  Electrons  are  heated  selectively  in  Ohmic 
heating,  and  since  the  transfer  of  electron  energy  to  the  ions  is 
inefficient,  this  technique  is  not  well  suited  for  the  heating  of  ions. 

(ii)  Collisional  Heating.  This  method  of  heating  is  often  referred 
to  as  magnetic  pumping.  The  magnetic  field  is  modulated  in  time, 

with  the  resulting  heating  in  the  perpendicular  direction  being 

transferred  through  collisions  to  the  parallel  direction.  The 

principle  involved  is  that  the  energy  increase  in  a  two  dimensional 

or  radial  adiabatic  compression  is  greater  than  in  a  three  dimensional 
3 

compression  .  This  technique  is  limited  to  low  frequencies,  because 
at  high  frequencies  only  the  surface  layers  are  compressed,  resulting 
in  only  the  surface  layers  being  heated. 

(iii)  Transit  Time  Heating.  In  this  method  of  ion  heating,  the 
magnetic  field  is  modulated  spatially  and  in  time.  This  process 
will  heat  the  ions  in  a  direction  parallel  to  the  magnetic  field, 
with  the  conversion  to  random  energy  resulting  through  collisions. 

The  basis  of  the  method  is  that  the  magnetic  moment  of  a  particle 
may  be  considered  to  be  constant  if  the  frequency  at  which  the 
magnetic  field  is  modulated  is  sufficiently  low. 

(iv)  Rapid  Compression.  Using  a  pulsed  magnetic  field,  the 
rapid  radial  compression  of  a  plasma  can  give  ion  energies  of 

400  (+)  ev.  However  the  conditions  required  for  a  positive  energy 
yield  in  a  thermonuclear  reaction  using  a  "fast-pinch"  technique 
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implies  an  instantaneous  release  of  approximately  10  joules  , 

which  is  equivalent  to  several  tons  of  TNT. 

(v)  ICRH.  A  technique  that  has  been  used  extensively  for  the 

heating  of  ions  in  a  plasma  is  ion  cyclotron  resonance  heating  (ICRH) . 

This  technique  has  been  used  for  ion  heating  in  the  B-65  stellerator^ , 

the  C  stellerator^ ,  the  B-66  mirror  device^,  and  in  various  other 

8  9 

mirror  machines  5  .  In  ICRH,  the  energy  in  an  ion  cyclotron  wave  is 
converted  into  the  random  motion  of  ions  through  the  process  of  cyclo¬ 
tron  damping.  Because  the  ions  are  accelerated  directly,  ICRH  schemes 
are  promising  for  the  attainment  of  high  ion  temperatures. 

In  the  schemes  that  have  been  employed  to  date,  an  ion  cyclotron 
wave  is  generated  in  a  region  of  the  plasma  where  is  typically  in 
the  range  of  .5  to  .85.  This  wave  is  then  allowed  to  propagate  along  a 
direction  of  decreasing  magnetic  field  strength,  until  the  local 
cyclotron  frequency  approaches  that  of  the  propagating  wave.  The 
wave  energy  will  be  absorbed  by  the  ions  in  this  region  through 
cyclotron  damping.  This  region  of  decreasing  magnetic  field  is 
referred  to  as  a  "magnetic  beach"  because  of  the  analogy  to  the 
breaking  of  water  waves  at  a  shoreline.  In  experiments  of  this  nature, 
the  electron  temperature  remains  much  less  than  the  ion  temperature 
because  of  the  effects  of  electron  radiation  and  the  selective  coup¬ 
ling  of  energy  to  the  ions  at  the  magnetic  beach. 

Three  techniques  that  have  been  used  for  coupling  the  rf  power 
from  an  external  source  to  an  ion  cyclotron  wave  in  the  plasma  are: 

(1)  Reverse-turn  induction  coil.  This  type  of  coil  is  normally 
referred  to  as  a  Stix  coil.  A  complete  analysis  of  its  use  in 
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coupling  rf  energy  from  an  external  source  to  a  plasma  is  given  by 

Stix^  (Chapter  5) ,  The  Stix  coil  has  been  used  for  the  generation 

of  ion  cyclotron  waves  in  the  B-65  and  C-slellerators ,  and  by  the 

11 

Kharkov  group  for  a  linear  discharge  in  a  uniform  field  ,  The 

efficiency  with  which  rf  energy  may  be  coupled  to  a  plasma  is  in  the 

12  3  13  3 

range  of  0,3  to  0„ 6  for  plasmas  with  densities  of  10  /cm  to  10  /cm  , 

Powers  up  to  1MW  at  25  M  Hertz  have  been  used^. 

An  example  of  the  results  that  may  be  obtained  for  the  heating 

12 

of  plasmas  by  the  use  of  a  Stix  coil  are  those  given  by  Nazarov 

13  3 

For  a  plasma  with  a  radius  of  3,5  cm,  and  a  density  of  10  /cm  ,  an 

ion  energy  of  2  Kev  was  obtained.  The  electrons  remained  much  cooler 

with  a  temperature  of  approximatley  30  ev„ 

The  reverse  coil  is  suited  to  the  generation  of  ion  cyclotron 

waves  in  plasmas  of  moderate  density  only,  that  is  for  densities  of 
1 A  3 

less  than  10  /cm  ,  It  cannot  be  used  for  the  efficient  generation 

of  the  short  wavelengths  required  for  high  density  plasmas.  That  Is, 

as  the  spacing  between  the  sections  of  the  coil  become  smaller  than 

the  coil  diameter,  the  electromagnetic  field  of  the  coil  decreases 

very  rapidly  as  one  moves  radially  inward  toward  the  plasma  column, 

2 

Also,  the  Kharkov  group  has  reported  that  the  use  of  a  Stix  coil 

at  power  levels  greater  than  10  to  20  KW  has  resulted  in  the  density 

of  the  plasma  rapidly  decaying  during  the  heating  pulse.  This  would 

appear  to  place  an  upper  bound  on  the  maximum  power  that  may  be  used, 

(2)  Coaxial  Electrodes,  Rf  power  may  be  applied  to  coaxial 

electrodes  at  the  end  of  a  mirror  device  for  the  generation  of  ion 

13 

cyclotron  waves  in  a  plasma  ,  The  theoretical  coupling  efficiency 
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may  be  estimated  by  considering  the  device  as  a  coaxial  waveguide, 

2  2 

filled  with  a  material  with  a  dielectric  constant  of  4ttc  /B  .  This 
approximation  is  valid  for  frequencies  much  lower  than  the  ion  cyclo¬ 
tron  frequency, 

14 

In  the  experiment  performed  by  Boley,  Wilcox  et  al  ,  the 

transfer  of  energy  from  a  1  MW,  8,3  M  hertz  oscillator  to  a  plasma 

12  3 

with  a  density  of  6  x  10  /cm  was  found  to  be  65%  efficient.  Better 
than  90%  of  the  wave  energy  was  absorbed  by  the  ions  at  the  magnetic 
beach. 

9  Iff 

A  similar  experiment  was  performed  by  Shvets  ’  on  the  Vikhr’ 

(Whirlwind)  device.  A  dense  plasma  was  produced  in  a  metal  discharge 

tube.  Ion  cyclotron  waves  were  then  generated  near  the  magnetic 

mirror,  and  propagated  in  the  axial  direction.  The  wave  energy  was 

absorbed  by  the  ions  through  cyclotron  damping  at  a  magnetic  beach. 

In  this  experiment,  an  ion  temperature  of  250  ev  was  obtained  in  a 

14  3 

plasma  of  density  of  2  x  10  /cm  .  Some  of  the  features  of  the  device 
9  15" 

used  by  Shvets  *  for  ion  heating  are: 

(i)  The  system  has  a  low  input  impedance  and  does  not  require 
high  voltages  for  the  introduction  of  rf  power  into  the  plasma. 

(ii)  Good  Coupling  between  the  plasma  and  the  electrodes  is  obtained, 

(iii) At  the  magnetic  beach,  the  energy  from  the  wave  is  initially 
coupled  to  ion  motion  which  is  perpendicular  to  the  static 
magnetic  field,  thus  improving  stability. 

(iv)  Two  gases  may  be  heated  simultaneously.  In  this  particular 
case,  argon  ions  reached  the  same  temperature  as  the  protons. 

The  use  of  a  coaxial  electrode  scheme  for  the  generation  of  ion 
cyclotron  waves  in  a  plasma  is  limited  to  mirror  confining  structures. 
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For  instance,  this  scheme  cannot  be  used  for  ion  heating  in  a  closed 
system. 

(3)  Spatially  Rotating  Magnetic  Field.  In  the  experiment  of  Karr, 

g 

Knapp,  and  Risenfeld  ,  a  plasma  was  injected  into  a  spatially  rota¬ 
ting  static  magnetic  field  from  a  coaxial  plasma  gun.  The  transverse 
component  of  the  magnetic  field  190  gauss)  was  produced  by  a 
helical  winding  placed  coaxially  in  the  mirror  device.  This  field 
was  superimposed  on  an  axial  field  of  approximately  3000  gauss.  By 
relating  the  coil  wavelength  (A^)  to  the  velocity  of  the  plasma  jet 
Oz)  by 

X  =  2ttv 
z  z 


the  plasma  saw,  in  its  frame  of  reference,  a  rotating  electromagnetic 

field  at  the  ion  cyclotron  frequency.  Through  the  ion  cyclotron 

damping  of  this  wave,  approximately  one  half  of  the  longitudinal  energy 

in  the  jet  was  converted  into  transverse  energy,  producing  a  plasma 

14  3 

with  T  =1  Kev  and  a  density  of  10  / cm  . 


1.2  Difference  Frequency  Harmonic  Ion  Heating 


In  an  attempt  to  overcome  some  of  the  problems  associated  with 
the  heating  schemes  discussed  in  Section  (1.1),  James  and  Thompson"^ 
have  suggested  the  use  of  the  difference  frequency  harmonic  which  is 
generated  by  the  nonlinear  mixing  of  two  waves,  for  ion  heating.  This 
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indirect  method  for  generating  low  frequency  waves  (oj  ^  in  a 
plasma  is  inherently  inefficient  because  the  waves  used  in  the  nonlin 
ear  mixing  process  are  selected  such  that  they  readily  penetrate  the 
plasma.  However,  since  the  difference  frequency  wave  is  generated  in 
ternally  in  the  plasma,  a  large  volume  of  the  plasma  may  be  heated 
simultaneously.  This  may  be  contrasted  with  the  absorption  of  the 
wave  energy  by  the  ions  at  a  magnetic  beach  for  the  case  where  either 
a  Stix  coil  or  coaxial  electrodes  are  used  to  couple  rf  energy  to  a 
plasma.  By  choosing  the  frequencies  for  the  waves  used  in  the  non¬ 
linear  mixing  process  to  be  greater  than  the  electron  plasma  and 
electron  cyclotron  frequencies  (see  Section  (1.2.1)),  ions  may  be 
heated  in  plasmas  confined  by  either  open  or  closed  systems.  This 
may  be  contrasted  with  the  use  of  coaxial  electrodes  or  a  spatially 
rotating  magnetic  field  for  ion  heating,  in  which  case  a  mirror 
device  (open  system)  is  required. 

In  this  thesis,  the  use  of  Whistlers  in  the  nonlinear  mixing 

process  will  be  investigated.  It  will  be  shown  that  plasmas  with 

12  16  3 

densities  in  the  range  10  to  10  /cm  may  be  heated  through  the 
nonlinear  mixing  of  whistler  modes  below  the  infrared  region. 

As  mentioned  above,  the  waves  used  in  the  nonlinear  mixing 
process  are  selected  such  that  they  readily  penetrate  the  plasma. 

Two  possible  choices  for  these  waves  are  listed  below. 

(i)  High  Frequency  Waves.  These  are  waves  with  a  frequency 

greater  than  or  equal  to  the  electron  plasma  and  electron 
cyclotron  frequencies.  This  category  includes  the  use  of 
lasers  in  the  mixing  process. 


' 
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(ii)  Whistler  Waves.  These  are  waves  with  a  frequency  in  the 
range  between  the  ion  and  electron  cyclotron  frequencies , 
and  which  propagate  approximately  parallel  to  the  static 
magnetic  field. 

The  power  absorbed  by  the  ions  from  the  difference  frequency 
wave  will  be  shown  to  be  proportional  to  the  fourth  power  of  the 
driven  electric  field  intensity.  Therefore,  the  successful  implemen¬ 
tation  of  difference  frequency  harmonic  ion  heating  will  require 
the  use  of  driving  waves  with  high  electric  field  intensities. 

The  results  for  the  nonlinear  mixing  of  high  frequency  waves 
and  that  of  whistlers  are  summarized  in  Sections  (1.2.1)  and  (1.2.2) 
respectively. 


1.2.1  H  F  Heating 


17  18 

Two  schemes  have  been  suggested  recently  ’  for  the  use  of  waves 
with  frequencies  greater  than  or  equal  to  the  electron  plasma  and 
cyclotron  frequencies  in  the  mixing  process.  In  both  schemes,  the 
energy  absorbed  by  the  ions  through  the  collisional  damping  of  the 
mixed  wave  is  optimized  by  a  resonance  in  the  magnitude  of  the  second 
order  fields.  This  is  achieved  by  allowing  the  mixed  wave  to  approach 
a  natural  mode  in  the  plasma.  In  the  above  cases,  this  was  chosen 
to  be  an  extra-ordinary  Alfven  wave  propagating  perpendicularly 
to  the  static  magnetic  field.  The  requirements  for  a  field  resonance 
place  stringent  restrictions  on  the  allowable  f luctuatuions  in  the 


. 
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angle  of  incidence  for  the  sources,  the  density  of  the  plasma,  the  static 
magnetic  field,  and  the  frequency  of  the  waves. 


1.2.2  Whistler  Heating 


The  nonlinear  mixing  of  whistlers  for  the  heating  of  ions  in  a 
magnetized  plasma  will  be  investigated  in  this  thesis.  In  order  that 
collisionless  damping  effects  be  included,  a  kinetic  analysis  will 
be  used  to  describe  the  heating  and  mixing  process. 

The  use  of  a  second  order  ion  current  resonance,  obtained  by 

allowing  the  difference  frequency  wave  to  approach  the  ion  cyclotron 

frequency  for  optimizing  the  power  absorbed  by  the  ions  from  the 

difference  frequency  wave  will  be  investigated.  It  will  be  shown  that 

this  technique  may  be  used  to  gain  an  order  of  magnitude  in  the  power 

17 

absorbed  by  the  ions  over  that  obtained  by  James  and  Thompson  ,  as 
well  as  the  relaxation  of  the  restrictions  on  the  allowable  density 
fluctuations  and  perturbations  in  the  angle  of  incidence  for  the  incident 
waves.  A  gain  of  two  orders  of  magnitude  is  obtained  in  the  sensitivity 
to  frequency  disturbances,  and  one  order  of  magnitude  in  the  sensitivity 
to  static  magnetic  field  fluctuations. 

The  use  of  the  second  order  ion  current  resonance  coupled  with 
the  field  resonance  effect  described  in  Section  (1.2.1)  will  also  be 
investigated.  It  will  be  shown  that  this  technique  may  be  used  to 
gain  two  orders  of  magnitude  in  the  power  absorbed  by  the  ions  over 
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that  obtained  by  James  and  Thompson*- ^ .  A  gain  of  two  orders  of 
magnitude  is  obtained  in  the  sensitivity  to  frequency,  density  and 
d-c  magnetic  field  fluctuations.  The  heating  process  is  less 
sensitive  by  two  orders  of  magnitude  to  perturbations  in  the  angle  of 
incidence  for  the  incident  waves. 


CHAPTER  2  NONLINEAR  INTERACTION  OF  WAVES  IN  A  PLASMA 


In  this  chapter,  a  kinetic  analysis  will  be  used  to  describe 
the  nonlinear  mixing  of  two  waves.  The  waves  will  be  assumed  to  have 
frequencies  in  the  whistler  frequency  domain.  Ion  and  electron  con¬ 
tributions  to  the  net  current  and  charge  perturbations  will  be  con¬ 
sidered.  The  selection  of  the  static  magnetic  field  will  be  such  that 
the  collisionless  damping  effects  on  the  incident  waves  are  negligible. 
A  cold  plasma  analysis  may  be  used  to  show  that  terms  of  order  of  the 
electron  collision  frequency  divided  by  the  electron  gyro-frequency 

(v  ,/Q,  )  are  introduced  by  the  collisional  effects  on  the  incident 
ei  - 

waves,  provided  the  incident  frequencies  are  not  too  close  to  a  gyro 
resonant  frequency.  For  typical  plasmas  under  consideration,  this  term 
is  of  the  order  of  10  In  the  following  analysis,  the  collisional 

damping  effects  on  the  incident  waves  will  be  neglected.  Collisional 
effects  will  be  included  in  the  second  order  induced  currents. 


2.1  Kinetic  Equations 


The  kinetic  approach  used  to  describe  the  mixing  and  heating 
process  begins  with  the  collisionless  Boltzmann  equation  expressed 
in  Lagrangian  coordinates.  That  is,  the  change  in  the  distrib¬ 
ution  function  with  time  is  written  in  a  coordinate  frame 
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moving  with  the.  zero’ th  order  trajectory  of  a  particle.  To  solve  the 
resulting  kinetic  equation,  the  method  of  orbit  integration  given  in 
Stix"^,  will  be  extended  to  include  second  order  fields  and  currents 
(see  Appendix  A),  Later  in  Section  (2,4),  the  corrections  accounting 
for  effects  of  collisions  between  the  species  will  be  discussed. 

The  Boltzmann  equation  for  the  j ’ th  species  may  be  expressed 
as  follows: 


+ 


Zj (  E  +  V*_§\  ■  ?£_J  - 
~b  v  Mj  C  y 


0 


(2.1) 


where  Z^e  =  magnitude  of  the  charge  on  species  j 
=  sign  of  the  charge  on  species  j 
m,  =  mass  of  species  j 

The  change  in  the  distribution  function  along  the  zero' th  order 
trajectory  of  a  particle  is  given  by 


(--) 

\dt  /o 


2>rJ 

dr 

.  dx 

"dYL 

di 

IX 

df 

+  v  •  ^ 

eCj  v/x 

B0  ■ 

,  ^ 

C 

(2 o  2) 


The  solution  to  equation  (2,2)  will  be  expanded  as  follows: 


. 
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r  +  - 
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i  (  Kf  V  -  OJ^b) 

■Pe“ 

1  (  K  ^  1  ¥  "■  ^ 

1  (  K  -Y  -  to  t) 

e  + 
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E 

0 

—  3 

_  B  . 

_  B0 

3± 

§3 

—  — 

J 

,,,(2.3) 

+  Conjugate  terms  +  higher  order  terms 
where  K  =  and  to  =  co^  - 

The  equilibrium  distribution  function  for  the  ions  and  electrons  is 

”1“  -  — 

represented  by  f^  and  f  respectively 0  The  first  and  second  order 
perturbations  to  the  distribution  function  are  given  by  f  _  and 

±  y  Z. 

-f“ 

f~  respectively. 

The  zero' th  order  trajectory  for  the  ions  and  electrons  is  given 
by 

d\v'  -  £ v/  x  x 

dh 


where  v'  =  velocity  vector  at  time  t' 


or  ,  depending  upon  the  species  being  considered. 


eB0 
Vi +c 


Q 


£■  Bo  I 

Wl_  cl 


e  =  sign  of  charge  on  species  being  considered. 

The  solution  to  the  above  equation  may  be  obtained  by  setting 
u~  =  v’  +  iev*  and  applying  the  boundary  condition  that  v’  =  v  for 
This  gives 


t  '  =  t  o 
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Vx  =  Vx  COS  (  Qt  -  S2t0  -  <?  Vy  s.VW«t~S2tO 

Vy  =  f  \/x  si'n  (  SU- S7*0  t  Vu  cos  £  S24-S2*')  /0  /N 

O  30  \  Z  O  4J 

=  Vz 

The  velocity  components  in  equation  ( 2 0 4)  may  be  integrated, 
and  with  the  use  of  the  boundary  condition  that  _r'  =  _r  for  t*  =  t, 

X1  -  -  Vx_  Si  Sit'  +  SVy  (l  -  COS  Sit)  +  X 

si  Tl 

y  1  -  -  €  Vx  (  1  -  COS  si  t)  -  Vl^_  Sm  SI  r  +  Cj 

^  52 

z'  =  -Vzr  +  z  o » o  (2o 5) 

r  -  £--t' 


The  collisionless  Boltzmann  equation  given  by  equation  (2»1), 
together  with  equations  (202)  and  ( 2 0  3 )  may  be  used  to  express  the 
change  in  the  equilibrium  distribution  function  along  the  zero' th 
order  trajectory  of  a  particle  as  follows: 

(d  --  O 

Mfc  'o 

The  most  general  solution  to  the  above  equation  is  of  the  type 


(2.6) 


where  _L  and  z  refer  to  directions  perpendicular  and  parallel,  to  the 

2  2  2 

static  magnetic  field,  and  v^  =  v^  +  v^o 
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For  the  specific  case  considered  in  this  thesis,,  f  will  be  taken 

o 

to  be  Maxwellian,  with  equal  temperatures  in  the  directions  perpendicular 
and  parallel  to  the  static  magnetic  field 0  That  is, 


exp  /  - 

^  £KTt 


,  ,  3  . 

(  v±  +•  V 


(2.7) 


The  velocity  gradient  of  f~  may  be  expressed  as  follows: 


>£L 

1)  V* 


1  VtJ 


± 


where 


(2.8) 


■t  +  r; 

By  using  equation  (2.4),  it  may  be  shown  that  f~,  fQ^,  and  f^_  will  have 

the  same  form  in  primed  and  unprimed  coordinates.  This  is  the  result 
2 

of  (v,[)  and  v’  being  constants  of  the  zero  order  motion. 


2.2  Solution  for  the  Incident  Fields 


Through  the  use  of  equation  (2.1)  to  (2.3),  the  equations 
governing  the  first  order  perturbations  induced  in  the  plasma  by 
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the  incident  fields,  may  be  expressed  as 


-  (e,  +  y  y 

cT” 


plus  a  similar  equation  for  f^ 

Through  the  use  of  the  following  Maxwell’s  equation 


V  x  E 


the  integration  of  equation  (2.9)  along  the  zero  order  trajectory  gives 


. . . (2.10) 


The  first  order  current  in  the  plasma  may  be  obtained  by  taking  the 

+ 

velocity  moments  of  v  f^.  The  dispersion  relation  for  the  first  incident 
wave  may  then  be  obtained  by  substituting  the  resulting  expression  for 
the  first  order  current  into  Maxwell's  equations.  A  similar  analysis 
is  applicable  for  the  second  incident  wave. 
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2  o  2  a  1  Dispersion  Relation  for  the  Incident  Waves 


The  results  given  in  Stix 
the  first  order  current:  as, 


10 


(Chapter  9)  may  be  used  to  express 


I,  -  vi 0  e  (  <  v^+>  -  >)  =  /f  <fK  tOpJ  ^  •  t;  ...  (2.11) 

K  7T 

■*irr 

t  ^  00  ± 

where  <1  V  -  \  >L  -P;  dv 

- 

i  =  192  and  denotes  the  first  and  second  incident  waves  respec¬ 
tively. 

n^  =  plasma  density/cm^ 

to  ,  =  plasma  frequency  for  species  k 

P  K. 

1c 

=  mobility  tensor  for  species  k 
=  sign  of  charge  on  species  k 

In  the.  following  analysis,  only  two  species  will  be  considered,  namely 
ions  and  electrons  (k  -  +,  -) . 

If  finite  L armor  radius  effects  are  neglected,  as  well  as 
collisionless  damping  and  finite  temperature  effects 
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Equations  (2.11)  and  (2.12)  may  be  used  to  express  the  total  first 
order  currents  as  follows: 


1  _ 


where 
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By  using  equation  (2.3)  and  Maxwell's  equations 


V  x  £f  = 

-1  )B 

c. 

V  x  B  = 

r  4- 
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c 

be  shown 
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Through  the  use  of  the  continuity  equation  and  Poisson's  aquation, 
the  above  equation  may  be  expressed  as  follows: 

-  (k:  •  Jr)  i<r] 

. o  o (2 . 15) 

where  i  =  1,2 

In  solving  equation  (2„15)  for  the  dispersion  relation,  it  will  be 
convenient  to  use  the  refractive  index  for  the  plasma,  which  is 
defined  by 

ml:  =  e.  k; 

,..(2, 16) 

If  the  incident  waves  are  assumed  to  be  Whistler  modes  with 
the  z-component  of  equation  (2.15)  gives 
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The  x  and  y  components  of  equation  (2.15)  may  now  be  expressed  in 
matrix  form  as  follows: 
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where 


L 

2 


fzi  r  ( 1  -  7 

«  -  i) -  ct.)  J 


A  non-trivial  solution  for  E 
This  gives 


i 
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and  E‘ 


implies  that  det 


]  =  0. 


-  2(1  -«f?)(l+ 

10( 


0 


. . . (2,19) 


The  solution  of  equation  (2,19)  will  give  the  dispersion  relation  for 
the  incident  waves. 


2.2.2  Reflections  at  the  Plasma-Vacuum  Interface 


The  assumed  geometry  of  the  plasma-vacuum  interface  is  given 
in  Figure  (2.1).  An  oblique  configuration  is  considered  so  that  the 
results  obtained  may  be  applicable  to  the  case  where  the  effects  of 
perturbations  on  the  slope  of  the  plasma-vacuum  interface  is  studied. 
A  similar  calculation  for  the  reflection  coefficients  has  been  made 
by  Unz^^  ,  Schmidt  ^  and  Sluijter^. 
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Figure  2,1  Geometry  of  Plasma-Vacuum  Interface, 

The  notation  used  in  Figure  (2.1)  is  as  follows: 
r 

=  propagator  for  the  reflected  wave 

s 

=  propagator  for  the  incident  wave 

=  propagator  for  the  transmitted  component  with  RH  elliptical 

polarization 

2t 

K  =  propagator  for  the  transmitted  component  with  LH  elliptical 
polarization 

i  =1,2  and  denotes  the  first  and  second  incident  waves  respectively 
The  notation  that  will  be  used  for  the  refractive  indicies  of  the 
plasma  to  the  waves  given  in  Figure  (2.1)  will  be  that  used  for  the 
propagation  vectors,  with  n  replacing  K.  For  example 
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lt 


The  transmitted  portion  of  the  wave  incident  on  a  plasma  vacuum 
interface,  as  depicted  in  Figure  (2.1)  will  decompose  into  two  com¬ 
ponents,  one  right-hand  and  the  other  left-hand  elliptically  polarized 
3 

(Rose  &  Clark  ) .  If  the  incident  frequency  is  taken  in  the  range 
between  the  ion  and  electron  gyro-frequencies,  the  LH  polarized  compon¬ 
ent  will  not  be  transmitted  in  the  plasma.  This  is  the  result  of  the 
strong  attenuation  of  the  wave  through  cyclotron  damping  near  the  ion 
cyclotron  frequency,  and  the  negative  value  for  the  square  of  the 

refractive  index  for  frequencies  greater  than  the  ion  cyclotron  freq- 
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uency  .  However,  the  RH  elliptically  polarized  component,  normally 
referred  to  as  a  Whistler  mode  may  be  made  to  penetrate  the  plasma 
by  a  suitable  choice  for  the  incident  frequency  and  the  static  magnetic 
field.  It  is  this  component  which  will  be  considered  in  the  nonlinear 
mixing  process. 

The  following  relations  may  be  obtained  from  Figure  (2.1): 
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The  solution  to  equation  (2.19)  may  then  be  written  as  follows: 
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j  =1,2  and  denotes  waves  with  a  RH  and  LH  elliptically  polarized 
electric  field  respectively . 

I  T  + 

lo  The  solution  for 


For  the  cases  that  will  be  considered s,  e 


jti 
i  ' 


it 

in  equation  (2„21)  will  be  obtained  through  an  iterative  procedure., 

j 

A  first  approximation  to  the  solution  for  rf:  is  obtained  by  setting 

.1. 

it  it 

ei  to  zero.  An  approximate  value  for  eh  ~  may  be  obtained  by  substitut- 

i  t 

ing  the  above  solution  for  rv:  into  the  expression  given  by  equation 

it  it 

(2.21)  for  en  o  A  better  approximation  to  the  solution  for  may  be 

obtained  by  substituting  the  above  value  for  into  the  equation  for 

)  o  The  accuracy  of  the  expressions  for  the  dispersion  relations 

2 

may  be  improved  by  further  iterations.  In  the  limit  that  (or /oj^)  <<1, 

it  2 

and  if  terms  of  order  (c^  )  l 2  are  neglected,  the  dispersion  relation 
for  the  incident  waves  in  the  plasma  may  be  written  as 
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where  the  value  for 


is  obtained  by  substituting 


e 

COS 


9\ 


+  (3-  2  ,)  of  ") 

u);  / 


,  into  the  expression 


given  by  equation  (2.21)  for  e^. 

The  required  expressions  for  the  propagation  vectors  may  be 
obtained  by  substituting  the  results  from  equation  (2.21)  (or  from 
equation  (2.22)  in  the  limit  that  (e^t)^/2  <<  1)  into  equation  (2.20). 

The  transmitted  field  components  will  now  be  expressed  in  terms 
of  and  E^  by  the  use  of  equation  (2.17)  and  (2.18). 
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where  1  =  1,2  and  denotes  the  first  and  second  incident  waves 
respectively, 

j  =1,2  and  denotes  waves  with  RH  and  LH  elliptically  polarized 
electric  fields  respectively. 
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Through  the  use  of  Maxwell’s  equations,  the  magnetic  field  components 
associated  with  the  incident  waves  may  also  be  expressed  in  terms  of 
E1C  and  E2t. 

IX  IX 
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where  i  =  1,2  and  denotes  the  first  and  second  incident  waves  res¬ 
pectively  . 

j  =1,2  and  denotes  waves  with  RH  and  LH  eliiptically  polarized 


electric  fields 

respectively. 

The  following  relations  may  be  obtained  by  the 

and  (2 . 24) : 
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The  incident  waves  may  be  taken  as  being  transverse  in  free  space, 
that  is 


!  (E*  •  O  - 
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. . . (2.26) 


where  iC?  =  u^/c  (sin  (G  +  9^)  ,  0,  cos  (9  +  9^)) 


=  propagation  vector  for  the  i’th  incident  wave  in  free  space, 
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i  =  1,2 

By  using  equations  (2.24)  and  (2.26),  the  following  relations  may  be 
derived  for  the  magnetic  field  components  in  free  space: 

s  Sv 

B,x  =  -  cos  (  B  +  &;  ) 

6,'cj  "  ^  /co%  (  &  +  )  fc";x 

Biz.  =  Sm  <  ©  +  e*  )  ...(2.27) 

s  s 

Both  E.  and  E„  are  known,  since  these  are  the  chosen  components  for 
ix  ly 

the  incident  waves.  The  electric  field  components  perpendicular  to  the 

g 

y-axis  for  the  incident  waves  in  free  space  will  be  denoted  by  E_^  . 

That  is , 


—  ^  — ^>/s 

-i-t  =  r*  +  tiz-  ...(2.28) 

By  taking  the  reflected  wave  as  being  transverse  and  using 
Maxwell’s  equations  to  express  the  relationship  between  the  magnetic 
and  electric  field  components,  the  following  relations  may  be  obtained 
for  the  reflected  wave. 
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(2  o  29) 


A  definition  similar  to  that  made  in  equation  ( 2  0  28)  for  the 

incident  waves  is  now  made  for  the  waves  reflected  from  the  plasma- 

r 

vacuum  interface 0  That  is  E,  is  defined  as  follows: 
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The  application  of  the  tangential  boundary  conditions,  namely, 
the  continuity  of  the  electric  and  magnetic  field  components  tangen¬ 
tial  to  the  plasma-vacuum  interface  yields 
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(continued  on  page  29) 
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The  transmitted  field  components  may  be  uncoupled  from  the  reflected 
field  components  by  combining  the  first  and  last  equations  and  the 
second  and  third  equations  obtained  by  the  expansion  of  equation  (2,31) 
By  defining 

a-  =  co  S  (  9^  -  d)  Cos  &  -  SI  (  ©?  -  ©) 
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the  equations  governing  the  transmitted  field  components  may  be 
expressed  as  follows: 
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The  solution  to  the  above  equation  may  be  readily  obtained  through 
the  use  of  Cramer's  Rule0 

(  M  i5-  kJ  1 4  i  fc  Kj  i2  t  ^ 

(  A  -  K>jg  KJ  i3  Eil)  o  o  o  (2  o  34) 
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where 


Nlil  N  ,'zj  -  nJ  ,-£  M/ 


The  reflected  field  components  may  be  found  by  substituting  the 
results  obtained  in  equation  (2^34)  into  equation  (2°31) „  This  gives 
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With  the  help  of  equation  (2o3),  the  percentage  of  the  power  in  an 
incident  wave  that  is  reflected  at  the  plasma-vacuum  interface  may  be 
shown  to  be  equal  to  100 

The  normal  configuration  that  will  be  assumed  for  the  incident 
waves  and  the  plasma-vacuum  interface  is  depicted  in  Figure  (202)„  The 
above  relation  for  the  reflected  power  will  be  used  when  examples  are 


considered  in  Chapter  3o 


' 
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FIGURE  2.2.  Wave  Vectors  and  Plasma- Vacuum  Interface.  Only  the  waves 
considered  in  the  mixing  process  are  shown  in  this  Figure. 
The  two  strongly  attenuated  LH  polarized  modes  have  not 
been  shown. 


A  resonance  in  the  magnitude  of  the  fields  associated  with  the 
mixed  wave  may  be  realized  if  the  propagation  vector  of  the  mixed  wave 
approaches  that  of  a  natural  mode  in  the  plasma.  In  this  particular 
case,  if  the  mixed  wave  is  taken  as  propagating  approximately  perpend¬ 
icular  to  B^,  the  natural  mode  is  an  extra-ordinary  Alfven  wave.  If  a 
field  resonance  is  to  occur  in  the  plasma  for  the  case  that  K^lliio’  t^ie 
required  angle  of  incidence  for  the  first  incident  wave  may  be  approxi- 
mated  as  follows  by  the  use  of  the  results  given  in  Spitzer  for  the 
dispersion  relation  for  an  extra-ordinary  Alfven  wave: 
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where 


(Ap_  (  1  +■  M-  /yVt+) 


The  propagation  vector  for  the  difference  frequency  wave  (See  Figure 


-32- 


(2o2))  is  given  by 

it  l-fc 

K  =  K >  -  Ks 

x  ...(2.37) 

The  object  of  this  thesis  is  to  investigate  the  use  of  this  difference 

frequency  wave  for  the  heating  of  ions  in  a  plasma. 


2.3  Second  Order  Fields  And  Currents 


In  this  section,  the  second  order  fields  and  currents  will  be 
obtained.  Through  the  use  of  equation  (2.1),  (2.2),  and  (2.3),  the 
equation  governing  the  second  order  perturbations  may  be  expressed  as 
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(2.38) 


The  velocity  moments  of  f~  will  give  rise  to  three  terms.  The 
first  term  will  be  an  induced  component  from  which  the  ion  energy 
absorption  will  be  calculated.  The  remaining  two  terms  correspond  to 
driving  terms.  The  contributions  to  the  second  order  distribution 
function  will  be  denoted  as  follows: 


+ 


The  first 


second  and  third  terms  on  the  RHS  of  equation 

+ 


(2.39) 


(2.39)  will 


denote  the  resulting  contributions  to  f ~  from  the  first,  second,  and 
third  terms  respectively,  on  the  RHS  of  equation  (2.38). 

The  procedure  that  will  be  used  to  solve  for  the  various  com- 


-33- 


+  + 

ponents  of  f  given  in  equation  (3„39)  is  outlined  below  for  f  „ 

3  3  Z 

Specific  reference  to  a  particular  species  will  not  be  made  since 
necessary  signs  and  factors  are  preserved  in  terms  as  e,  Z,  and  m„ 

From  equation  (2.38) 

^3e  =  -  Ze_C  C  (  Fg  +  y  x  d±! 

T-*  00  )  ~c.  J  =..(2.40) 

-T 

By  using  the  expression  given  for  f ^  in  equation  (2.10),  the  above 
equation  may  be  expressed  in  the  following  form: 

=  *)  C  dr  dr'  G  (  ET*  (r-,r)  rO,  v)  .0.(2, 41) 

r  V' 

where 

with  the  integration  limits  on  t  and  x*  being  [0,°°)  and  [ 0 , 00 ] 
respectively.  The  derivation  of  equation  (2„41)  is  given  in  Appendix  B. 


2.3.1  Velocity  Moments  of  f^ 


Following  the  procedure  outlined  in  equation  (2.38)  to  (2.41) 
and  Appendix  B,  the  first  step  in  solving  for  f^?  is  to  find  an  expres¬ 
sion  for  f1(t’).  This  may  be  obtained  from  equation  (2.10).  By  using 
equation  (2.5),  the  resulting  phase  term,  i (K^ »_r '  ’  -  to^t*’)  may  be 
expressed  in  terms  of  the  variables  at  r’  and  t*. 


That  is 


. 
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i  Kj^  •  T  "  -  i  u),  t '  =  i  •  Y '  -  i  Ldx  t'  +-  |'4  (  -  l<t^  smlZr') 
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. . . (2.42) 


where  x ’  =  t ’  -  t ’ ' 

By  performing  a  similar  analysis  for  the  term 
i  £  \  .  y*  — 


£i  e 
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f (_r '  , v '  , t 1  )  may  be  expressed  as  follows: 
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From  the  results  given  by  equation  (2.8)  for  a  plasma  with  a  Maxwellian 

equilibrium  distribution,  and  with  equal  temperatures  in  the  directions 

perpendicular  and  parallel  to  the  static  magnetic  field,  it  is  easy  to 

show  that  (f  -  v'  f  )  =  0.  For  this  particular  case,  the  above 
oz  z  o 

expression  for  f  ^ (_r  ’  ,v;’ ,  t ' )  simplifies  to 
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..,(2. 43) 

By  using  that  =  c/w^K^  x  ,  equation  (2*,  40)  may  be  expressed 
as  follows: 
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where  f^(r’  ,v'  ,  t')  is  given  by  equation  (2043), 

If  only  the  transmitted  incident  waves  with  a  RH  elliptically 
polarized  electric  field  are  considered  in  the  mixing  process,  the 
following  relations  may  be  obtained  from  equation  ( 2  0  23 ) : 


*  * 
Eey  =  -  • 


..,(2, 45) 


In  deriving  equation  (2„45),  the  propagation  vector  for  the  second 
incident  wave  was  assumed  to  be  parallel  to  the  static  magnetic 
field  in  the  plasma  (see  Figure  2o2)„  Through  the  use  of  equations 
(2„43)  and  (2,45),  the  partial  derivative  of  f  (r' ,v? , t ' )  with  respect 
to  vT  may  be  shown  to  be  equal  to 
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where  v!  =  v'  +  iv* 

+  x  y 

By  the  use  of  equations  (2.4)  and  (2.5),  the  variables  in  equations 
(2.44)  and  (2.46)  may  be  expressed  in  terms  of  _r,  v,  and  t,  where 
T  =  t-  t '  .  Equations  (2.44)  and  (2.46)  may  then  be  combined  to 
give  the  following  equation  for  f  : 
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B^q  ,  see  Figure  (2„2)) 


=  difference  frequency 
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The  next  step  in  the  analysis  will  be  to  express  the  integrand  in 
equation  (2.47)  in  terms  of  the  various  powers  of  the  velocity  compon¬ 
ents  v^  and  Vo  Equation  (2.47)  may  then  be  expressed  as: 
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a  -  -  k*  sivt  C  r  +  x' ) 
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b  =  (  1  -  CO<L  n  LX  +  r' )) 
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^  Kzvz)  T  +  (cJL  -  k^v?)r' 


In  the  next  step,  an  integration  over  the  velocity  components 
v^  and  will  be  performed  for  the  various  velocity  moments  of  f^0 
Following  the  notation  used  in  Stix“^,  the  velocity  moment  for  q(v)  in 
the  x-y  velocity  space  is: 
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..,(2. 50) 


Note  that  the  expression  for  f  ^  given  by  equation  (2=48)  still  requires 
an  integration  over  t  and  x'.  This  will  be  performed  later  in  the 
analysis o 

In  the  integration  of  equation  (2*48)  over  v^  and  v  ,  integrals 
of  the  following  type  appear  (see  Stix^^,  pp  173): 
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By  completing  the  square  in  the  exponential, 
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A  standard  integral  form  is  obtained  by  defining  u  =  v  -  i  akT  „  The 
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results  for  p  =  0,  1,  2,  and  3  are  as  follows: 


(  Dj  a,) 

<5  (  1,  a) 
G  (  2,  a) 
(3  (  3,  a) 


. . .  (2.52) 


When  an  integration  over  and  v  is  performed,  it  is  evident  from 

2  2 

equation  (2.52)  that  an  exponential  term  with  the  exponent  -(a  +  b  )KT/2m 

will  appear  for  all  values  of  p.  This  exponent  may  be  expanded  by  the 
use  of  the  definitions  for  a  and  b  to  give 


-  (  a2  4-  b  )  KJ  =  -  (  1  -  Cos RC'i'+  r’) ) 


. . . (2.53) 


where  "A  - 

Qe  y« 

The  factor  A  is  proportional  to  the  square  of  the  ratio  of  the  Larmor 
radius  to  the  wavelength  perpendicular  to  the  static  magnetic  field. 
This  factor  represents  the  finite  Larmor  radius  correction  terms  and 
will  later  appear  as  an  expansion  parameter.  In  the  case  being  con¬ 
sidered,  namely  an  electromagnetic  wave  with  a  frequency  somewhere  in 
the  range  between  the  ion  and  electron  cyclotron  frequencies,  incident 
on  a  plasma  half-space,  (See  Figure  (2.2))  the  following  inequality 
may  be  shown  for  A  : 
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. . . (2.54) 
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where  T  =  ion  temperature  in  degrees  Kelvin 


COi  =  Cx 


*/a 


Q-h  =  (S7.4<T2._)  =  Hybrid  frequency 


If  the  ions  and  electrons  are  taken  to  be  at  the  same  temperature, 

-r  -  («±) 

(  k)+/ 

By  using  the  notation  given  in  equation  (2.62),  the  various 
velocity  moments  of  f  in  the  x-y  velocity  space  may  be  expressed 
as  follows: 
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^  A  COSfl  ( T  +  t’)  .  1  T  ,  1  ^ 

The  term  e  will  be  expanded  by  the  use  of  the  Bessel 

function  equality 
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where  I  (A)  =  i  n  J  ( i A) .  In  order  that  terms  with  a  common 
n  n 

exponential  dependence  may  be  collected  so  that  the  integrations  over 

x,  x ' ,  and  v  may  be  performed,  the  coefficients  appearing  in  equation 
z 

2 

(2.68),  namely  a,  a  ,  ...  must  be  evaluated.  This  is  done  below. 
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where  I  ^4  kt 

r"  =  r  r  r' 

In  the  next  section,  the  integration  with  respect  to  i'  in  equation 
(2.55)  will  be  performed.  This  requires  that  the  terms  with  a  common 
exponential  dependence  in  t'  be  collected.  After  some  manipulation, 
equation  (2.55)  may  be  expressed  as  follows: 
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The  integration  with  respect  to  x'  in  equations  (2  ^  58)  to  ( 2 0 60) 
will  now  be  performecL  The  asymptotic  form  of  the  assumed  solutions 
for  the  incident  fields  is  as  follows: 
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However,  in  the  limit  that  t'  ■+  °°,  the  field  must  vanish  since  it  must 
be  switched  on  at  some  finite  time,,  Therefore,  co^  may  be  thought  of  as 
having  a  very  small,  but  finite  positive  imaginary  component .  Solutions 
for  steady  state  and  damped  oscillations  may  be  obtained  from  the  sloutions 
for  which  Im(aj^)  >  0,  by  an  analytic  continuation  in  the  frequency 
domain  (see  Stix^,  pages  169  and  179) .  Two  types  of  integrals  involving 
x’  occur  in  equations  (2„58),  (2„59)  and  (2*60). 
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(i)  Integrals  with  x’  in  an  exponential  term  only.  These  are 
integrals  of  the  type 

OO 
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Jo 

If  10^  has  a  small  positive  imaginary  component,  this  integral 
may  be  shown  to  be  equal  to 

i  C  )/  (  ^  SI  +  K 1  z. 


(ii)  Integrals  involving  x’  in  an  exponential  term  and  linearly  in 
the  remainder  of  the  integrand.  These  are  integrals  of  the 
type 
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o 

The  above  integral  may  easily  be  integrated  by  parts,  and  if 
is  taken  to  have  a  positive  imaginary  component,  the 
application  of  L'Hospital's  rule  will  give  the  following  result 
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The  results  obtained  in  parts  (i)  and  (ii)  may  be  expanded  asymptoti¬ 
cally  as  follows: 
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The  expansions  given  by  equation  (2,61)  will  converge  if  co^  is 

sufficiently  removed  from  the  ion  and  electron  cyclotron  frequencies 

and  their  harmonics .  The  integration  with  respect  to  x’  in  equations 

(2.58)  to  (2.60)  will  now  be  performed  and  terms  with  the  exponent 

2 

i(u)  +  nfi  -K  v  )  summed.  Terms  of  the  order  (v„/v  .  )  will  be 
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neglected,  where  v  is  the  thermal  velocity  for  the  species  being 

considered,  and  v  ,  and  v  _  are  the  phase  velocities  along  the  static 

plz  p2 

magnetic  field  for  the  first  and  second  incident  waves  respectively  in 
the  plasma.  Through  the  use  of  the  following  definitions: 
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equations  (2.58)  to  (2.60)  may  be  rewritten  as  follows  after  the  inte¬ 
gration  with  respect  to  x1  is  performed: 
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The  integrals  still  to  be  evaluated  are  of  the  type 

a C  c^> 

Fp  =  K^  C  d  vz  f  dT  Vj,  exp  f  -  kW^  +  ('2  (+*>  +tjX2  -  K^V^))  .  .  .  (2 . 66) 
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-  jt>  o 

where  p  =  0,  1,  2,  or  3  .  The  results  given  in  Stix^  for  integrals 
of  this  type  may  be  generalized  to  include  the  case  where  p  =  3.  In 
order  to  demonstrate  the  convergence  of  the  integration  with  respect 
to  x  in  equation  (2.66)  for  either  positive  or  negative  values  for  the 


-49- 


imaginary  component  of  co,  the  square  in  the  exponent  may  be  completed 
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As  is  evident  from  equation  (2.67),  convergence  will  exist  for  either 
positive  or  negative  values  for  Im(w).  Therefore,  F^  may  be  considered 
to  be  an  entire  function  of  w. 

In  the  next  step,  F^,  F^,  and  F^  will  be  expressed  in  terms  of  Fq» 
For  Irn(m)  >  0,  the  integration  with  respect  to  r  in  equation  (2.66) 
may  be  performed  first,  to  give 

F?  =  I  C  dVz  \l ^  e.  g‘,<'r  .  ..(2.68) 
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Since  F  is  an  entire  function  of  w  (see  equation  (2„67)),  the  value 
P 

for  in  the  lower  w  half-plane  may  be  obtained  by  analytic  contin¬ 
uation. 

From  equation  (2.68) 
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The  integrand  in  the  above  equation  may  be  expressed  as  follows: 
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The  integration  with  respect  to  v  may  then  be  performed  to  give 
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In  a  similar  fashion 
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Either  an osymp tot ic  or  convergent  series  may  be  written  for  F  , 
the  latter  being  required  when  oo  is  close  to  a  multiple  of  the 
respective  gyro-frequency  or  when  the  Landau  damping  of  a  wave  is 
significant.  The  convergent  series  for  Fq  Is  obtained  by  first 
integrating  with  respect  to  v  In  the  expression  given  by  equation 

Za 

(2.67)  for  F  .  By  defining 
cU  n  SI )  [~FT 

K  ^  2KT  O  o  o  (2  O  71) 

and  making  the  following  transformation  of  variables, 
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the  Integration  with  respect  to  x  may  be  performed.  This  gives 
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where  S(z)  is  the  complex  error  function  and  is  given  by 


Slz) 


fe^dt 


. . . (2.73) 


The  function  F^  at  this  stage  may  be  evaluated  in  two  ways.  One  is  the 
use  of  tabulated  values  given  by  B.D.  Fried  and  S. D. Conte for  the 
function 


=  IT  ^  {rr  exp  -Ct2)  +  2i  S(z^  J 


Another  is  the  development  of  either  a  convergent  or  asymptotic  series, 
depending  upon  the  absolute  value  of 

A  convergent  series  is  obtained  by  the  integration  of  equation 
(2o73)  by  parts.  This  yields 
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This  procedure  may  be  repeated  to  give  the  series 
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Convergence  will  exist  for  \z\  <1,  which  implies  that  if  the  conver¬ 
gent  series  is  to  be  used,  lanl  < 

An  asymptotic  expansion  for  F  ,  valid  when  |a  |  >  1,  is  given  in 
Stix^  (pp.179).  The  technique  used  in  obtaining  this  expansion  will 
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be  briefly  outlined  below.  From  equation  (2,68) 
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where  the  symbol  P  denotes  the  principal  value  of  the  integral. 

If  the  expression  for  F  in  equation  (2,69)  is  compared  with  equation 

O 

(2,75),  an  alternate  expression  for  S(z)  is  obtained,  namely 


-L 

Six.)  =  _i_  pP  e  dt 

—  G<> 


,,,(2,76) 


The  asymptotic  expansion  for  S(z)  may  be  obtained  by  first  expanding 
the  integral  given  in  equation  (2,76)  for  z  real,  then  forming  the 
analytic  continuation  info  the  complex  plane  using  a  Taylor  series. 

The  asymptotic  expansion  for  S(z)  in  the  region  where  j Re (a  ) | > | Im(an) j 
is  given  by 
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In  the  region  where.  |  Re  (a^)  j  <  j  Im(a  )  |  ,  a.  different  expansion  is 

required.  Thus  the  asymptotic  expansion  for  5(z)  exhibits  a  Stokes 

10  25 

phenomenon  (see  Stix  (pp  180)  and  Morse  and  Feshbach  ) ,  An  asymptotic 
expansion  in  this  region  will  not  be  required  in  this  thesis. 

Equations  (2,63)  to  (2,65)  will  now  be  used  for  obtaining  expressions 
for  the  electron  and  ion  driving  current.  The  remaining  Integra- 
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tions  with  respect  to  x  and  v  in  these  equations  may  be  evaluated 

Z 

in  terms  of  F_ ,  F^,  F^,  and  F^.  F^  will  be  approximated  as  follows: 
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Where  is  given  by  equation  ( 2 o  71) .  In  the  second  expression  in 
equation  (2.74),  the  first  term  gives  the  contribution  of  the 
collisionless  damping  effects,  the  first  component  in  the  second  term 
gives  the  normal  cold  plasma  contribution,  while  the  second  component 
gives  the  first  temperature  correction  term. 

In  obtaining  the  various  current  components,  only  terms  to  the 
zero’th  order  in  X  (where  A  is  given  by  equation  (2.53))  will  be 
included.  Such  an  analysis  will  be  sufficiently  general  to  account 
for  Landau  or  Cherenkov  damping  and  cyclotron  damping  at  the  cyclotron 
frequency  for  the  species.  Expansions  to  higher  orders  in  X  may  be 
used  to  account  for  finite  Larmor  radius  effects  as  well  as  transit 
time  damping,  that  is,  the  effects  of  fluctuations  in  the  static 
magnetic  field.  Such  a  generalization  would  also  account  for  the 
effects  of  cyclotron  damping  at  the  higher  harmonics  of  the  gyro- 
frequency. 

From  equation  (2.39)  it  is  evident  that  the  second  order  current, 
given  by  the  velocity  moments  of  f^  will  have  three  separate  compon¬ 
ents,  two  driving  components  and  one  induced  component.  The  notation 
that  will  be  used  to  represent  the  driving  terms  is  as  follows: 


. 


-54- 


J 


t  (k) 

3a 


Ck) 


=  n. 


'eCt  ^3<)  =  A  ^  FT;  (n>| 


n.--  °£> 


(2.79) 


where  i  =  1,2,3  and  denotes  components  in  the  x,  y,  and  z  directions 
respectively. 

k  =  1,  2  and  complies  with  the  notation  used  in  equation  (2,39) 

n  refers  to  a  used  in  the  evaluation  of  F  ,  F..  ,  F„,  and  F„. 
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The  II  coefficients  are  evaluated  to  order  zero  in  A~  through 
the  use  of  equations  (2.63)  to  (2.65),  and  are  enumerated  below 
for  k  =  20  Only  the  non-zero  terms  are  listed. 
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2.3.2 


Velocity  Moments  of 


f 


31 


In  this  section,  the.  distribution  function  f 
moments  will  be  evaluated.  From  equations  (2.38) 
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^  and  its  velocity 
and  (2„ 39)  , 


...(2- 81) 


In  the  analysis  to  follow,  specific  reference  to  ions  or  electrons 
will  be  dropped  until  currents  are  to  be  evaluated.  The  first  step 
in  the  analysis  will  be  the  evaluation  of  f  (t')o  This  is  given  by 
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The  assumed  exponential  dependence  for  the  second  incident  wave  is 
as  follows: 

_ k  N  -!  (  Ka'l"  -  ^2.  4'0 

t  2.  ( !■"/  =  £”2  £ 

Through  the  use  of  the  Maxwell  equation 


Vx  t  -  -A 

c  ~dt 

equation  (2.82)  may  be  written  in  terms  of  alone.  With  the  aid 
of  equation  (2.4)  and  (2.5),  the  variables  in  _r 1  '  and  t'*  may  be 
expressed  in  terms  of  and  t’.  If  the  equilibrium  distribution 
function  (f  )  is  assumed  to  be  Maxwellian  with  equal  temperatures 
in  the  directions  perpendicular  and  parallel  to  the  static  magnetic 
field,  and  if  only  the  RH  circularly  polarized  portion  of  the  second 
incident  wave  is  considered  in  the  nonlinear  mixing  process 
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After  the  velocity  gradient  and  dot  product  in  equation  (2.81) 
are  performed,  the  variables  in  _r*  and  t’  may  be  expressed  in  terms 
of  r_  and  t  by  the  use  of  equations  (2„4)  and  (2.5).  The  justification 
for  this  step  is  discussed  in  Appendix  B.  If  the  propagation  vector 
is  taken  to  be  parallel  to  the  static  magnetic  field  (see  Figure 
(2.2)),  x1  will  appear  uncoupled  from  the  variable  x  in  the  exponent 
of  the  equation  for  f  .  This  may  be  compared  with  the  exponent  in 
equation  (2.47)  which  involves  terms  such  as  cos  ft  (x  +  x’).  This 
presents  a  considerable  simplification  over  the  analysis  for  f  , 
since  now  the  integration  with  respect  to  x ’  may  be  performed  before 
any  of  the  velocity  integrations  or  the  integration  over  x.  Two  types 
of  integrals  are  encountered  in  the  integration  with  respect  to  x ' . 
These  are 
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I£  =  C  d'ir'  (  )  t'  exp -i  (  <Pg*  _  KgVa  +  zo)  r '  ...(2.84) 


where  (  )  is  independent  of  x’.  The  integration  may  readily  be 

carried  out  by  parts.  If  is  assumed  to  have  a  positive  imaginary 
component,  integrals  1^  and  become 
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The  above  results  are  expanded  by  the  use  of  an  asymptotic  expansion 
as  in  equation  (2.61). 
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The  expansion  will  converge  if  the  frequency  u  is  sufficiently 

removed  from  the  ion  and  electron  cyclotron  frequencies  and  their 

harmonics.  As  in  the  analysis  for  the  velocity  moments  of  f  ,  terms 

2 

of  the  order  of  (v„  /v  .  .)  will  be  neglected,  where  vQ  is  the 

U  p  XZ  ^  Z ) 

thermal  velocity  for  the  species  being  considered,  and  v  ^  and  v  ^ 
are  the  phase  velocities  along  the  static  magnetic  field  for  the 
first  and  second  incident  waves  respectively,  in  the  plasma. 

The  integrand  in  the  equation  for  f  ^  will  now  be  expressed  in 
terms  of  the  various  powers  of  the  velocity  components  v^  and  v^  (see 
equation  (2.48)  for  a  similar  equation  for  f ^)  • 
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If  <J)  is  defined  as  follows: 


dj  =  -  "A  (  i  -  cos  S2t)  +  i  (  u>  -  Kzln/z')  ^ 


.0.(2. 88) 


the  results  from  equations  (2.51)  to  (2.53)  may  be  used  to  express  the 

velocity  moments  of  f0,  over  v  and  v  as  follows: 
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In  order  that  terms  with  a  common  exponent  can  be  collected  so  that 

the  remaining  integrations  over  v^  and  x  can  be  performed,  the  various 

2  2 

coefficients  such  as  a  (kT/m)  appearing  in  equations  (2.89)  to  (2.91) 


must  be  evaluated.  These  are  obtained  by  replacing  (x  +  x')  by  x  in 
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equation  (2.57).  By  the  use  of  the  Bessel  function  equality  given  by 
equation  (2.56),  the  velocity  moments  of  f^  given  by  equations  (2.89) 
to  (2.91)  may  be  expressed  in  an  infinite  series  as  follows: 
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The  integration  with  respect  to  t  and  v  ,  still  to  be  performed 
in  equations  (2.92)  to  (2.94)  are  of  the  same  type  considered  in  equation 
(2.66)  for  the  velocity  moments  of  f^*  By  using  the  results  obtained 
for  the  integral  in  equation  (2.66),  and  the  notation  for  the  driving 

currents  as  given  by  equation  (2.79),  the  II  coefficients,  evaluated 

+ 

to  order  zero  in  A  for  k  =  1  may  be  expressed  as  follows: 
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Equations  (2.79),  (2.80),  and  (2,95)  may  now  be  used  to  express 
the  total  second  order  driving  current  with  the  exponent  i(K._r  -  wt) 
as  follows: 
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where  k  =  1,  2  (see  equation  (2,39)),  If  the  ions  and  electrons  are 
assumed  to  be  at  the  same  temperature,  the  total  second  order  driving 
current  may  be  expressed  as  follows: 
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where  rr,,^  -  rr.,„i  +  rr.;^ 

i  =  1,  2,  3  and  denotes  components  in  the  x,  y,  and  z  directions 
respectively. 

+ 

The  coefficients  are  evaluated  below  through  the  use  of  equations 


(2.80)  and  (2.95). 
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2.3.3  Velocity  Moments  of  f 


The  results  for  the  induced  second  order  current,  namely  the 
current  that  is  obtained  from  the  velocity  moments  of  f  in  equation 
(2.39),  may  be  obtained  directly  from  the  results  given  by  Stix^  ( 
section  9.2).  In  the  analysis  by  Stix"*^,  the  currents  are  obtained  by 
first  defining  a  mobility  tensor  as  follows: 


' 


t 
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(2.99) 


10 


As  in  Stix  (section  10,8),  a  correction  to  the  mobility  tensor  will  be 

made  to  account  for  collisional  effects  between  the  species  under 

consideration.  Before  this  can  be  accomplished,  the  mobility  tensor 

for  a  collisionless  plasma  must  be  obtained.  Terms  of  order  A  will  be 

neglected  in  the  evaluation  of  the  mobility  tensor  below.  The  terms 

Fq,  Fr  ...  (where  is  given  by  equation  (2.66))  will  be  denoted  as 

F  /  N ,  F- ,  N ,  ...  where  n  refers  to  the  a  used  in  evaluating  these 
o(n)  l(n)  n 

terms.  The  various  terms  in  the  mobility  tensor  are  written  so  that 
they  will  be  applicable  to  ions  or  electrons  by  the  appropriate 
substitutions  for  e,  m,  T,  and  ft. 


(2.100) 
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where  is  the  tensor  component  relating  the  velocity  component  in 

the  i'th  direction  to  the  electric  field  component  in  the  j ’ th 
direction^  By  the  use  of  equations  (2,97)  and  (2,99),  the  total 


second 

order 

current  may 

^3  - 

a 

An'Tl 

where 

=  m;  -ds 

J  id 

=  in  hd 

+  _J . 
4rr 


.  (2.101) 


U) 


2.4  Correction  to  Mobility  Tensor  for  Ion-electron  Collisional  Effects 


The  mobility  tensor  will  now  be  corrected  for  ion-electron 

collisional  effects.  When  a  resonance  in  the  second  order  current 

is  obtained  by  allowing  the  difference  frequency  co  to  approach  the  ion 

cyclotron  frequency,  collisional  effects  will  be  relatively  unimportant^ 

if  v .  <<  fkT  K  (where  v.  =  m  /m,  v  „).  This  may  be  seen  by 

le  \ —  z  ie  -  +  ei 

v  m 


considering  the  mobility  tensor  component  .  For  a  collisionless 

jXX 

plasma , 
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The  expression  that  is  obtained  for  from  a  cold  plasma  analysis, 

JXX 

with  collisional  effects  included,  is 


(Cold')  -  iQ+ 
2 


fu- s iS) 


-  (cj+f2-l  (cJ-Sl-f)  +  iu>£  (u)-  Q-i)  (oJ+Qi)  4  iu>£  . 


. . . (2.103) 


where  v  =  v  _^(1  +  m_/m+) .  As  to  approaches  the  ion  cyclotron  frequency, 

a  resonance  in  both  the  collisionless  and  cold  plasma  expressions 

for  is  obtained.  However,  the  magnitude  of  this  resonance  will 

1/2 

be  bounded  by  cyclotron  damping  effects  for  (to  -  )  <  (2KT/m  )  K 

I  z 

in  the  case  of  a  collisionless  plasma  (equation  (2.102),  or  by 

collisional  effects  for  (to  -  fi+)  <  v  (equation  (2.103))  in  the  case 

of  a  plasma  where  collisional  effects  are  dominant.  Therefore,  if 

+  1/2 

v  <<  (2KT  /m+)  ,  cyclotron  damping  will  be  the  dominant  mechanism 

for  ion  energy  absorption  from  the  difference  frequency  wave. 

+  1/2 

For  the  case  that  v.  >>  (2KT  /m  )  K  ,  the  absorption  of  energy 

1 01  !  Z 

by  the  ions  from  the  second  order  fields  will  be  determined  exclusively 

by  collisions.  This  will  occur  for  very  dense,  low  temperature  plasmas. 

However,  for  this  case  there  is  no  need  for  a  kinetic  analysis,  since 

a  cold  plasma  analysis  will  suff ice26o  This  case  will  not  be  considered 

in  this  thesis,  since  the  use  of  cyclotron  damping  presents  a  much 

more  powerful  technique  for  the  heating  of  high  temperature  plasmas. 

Outside  the  "inner  resonance  region','  that  is  for  (to  -  fi+)  > 

+  1/2 

(2KT  /m+)  K^,  collisional  effects  will  be  included  in  the  mobility 

tensor  by  the  use  of  the  cold  plasma  theory  as  described  by  Stix"^ 

+  1/2 

(section  (10.8)).  That  is,  for  the  case  that  v.  <<  (2KT  /m  )  K  , 

1 0  i-  z 

collisional  effects  will  only  serve  to  broaden  the  frequency  range  for 


:  U 
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resonant  ion  energy  absorption*, 

In  the  cold  plasma  analysis  of  the  nonlinear  mixing  of  HF  waves 
17 

by  James  and  Thompson  ,  collisional  effects  introduce  terms  of  the 

order  <  v.  /  (oo  -  fl.)  into  the  expressions  for  the  second  order 

driving  currents*,  The  net  effect  of  these  terms  in  the  expression  for 

the  ion  power  absorption  from  the  second  order  fields  is  of  the  order 
2 

(v_^e/(w  -  fi+))  .  For  the  case  of  whistler  heating,  a  cold  plasma 

analysis  will  give  similar  results,  namely  that  collisional  effects 

introduce  terms  of  the  order  <  (v,  /  (oo  -  fl,))  into  the  second  order 

driving  currents .  To  check  if  similar  results  are  applicable  to 

whistler  heating,  the  cold  plasma  equations  were  programmed  on  the 

IBM  360  computer.  An  error  of  one  percent  in  the  ion  power  absorption 

was  introduced  when  collisional  effects  were  neglected  in  the  driving 

16  3  6  o 

terms  for  a  plasma  with  a  density  of  10  /cm  ,  a  temperature  of  10  K, 

and  with  (oo  -  ft  )  -  10  v .  .  This  is  what  the  results  for  the  nonlinear 

+  le 

mixing  of  HF  waves  predict. 

The  method  used  by  Stix^  for  including  collisional  effects  in 
the  mobility  tensor  is  outlined  below.  The  analysis  begins  with  the 
conservation  of  momentum  equation  for  species  j. 


n  ■  m  •  A  v/j 


+  Vj 

c 


) 


+ 


ZijK  '  (  V*  -  Xi) 

K 


, e . (2.104) 


where  the  last  term  on  the  RHS  of  the  above  equation  represents  the 
effect  of  collisions  with  species  k.  The  effect  of  collisions  will  be 
accounted  for  by  introducing  an  additional  velocity  component  Av^ .  By 
the  use  of  the  mobility  tensor  notation  given  by  equation  (2.99),  the 
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introduction  of  A_v^  into  equation  (2,104)  yields 


-i  ^  +  A  ^  ^  •  E  - 


from  which  AM^  may  easily  be  shown  to  be 


A 


hi  'A  1  ’  =jK ' 


K 


where  v.,  =  v„,  I  .  If  the  conservation  of  momentum  is  to  hold  for 
-jk  jk  = 

the  various  species,  the  following  relation  must  hold. 


ni  ^jK 


^yl  jr/7  A 

'K  K  Kj 


. . . (2.105) 


Eor  the  case  under  consideration,  namely  for  a  plasma  composed  of 
two  species  (ions  and  electrons) 


4  £4  3 


fie.  M  J  •_!  •  (  U3  -  H  j) 


A  Ms 


—  £  a.  I  M,  •  1  *  (  M 

S2+  ~ 


. . . (2.106) 


where  the  mobility  tensors  on  the  RHS  of  equation  (2,106)  are  those 
obtained  by  neglecting  collisions.  The  technique  used  above  may  be 
shown  to  be  equivalent  to  a  cold  plasma  analysis  in  which  just  the 


first  two  terms  in  the  asymptotic  expansions  of  _ 1 


f  1  +■  (  1+  Yrt./y/l-t-')  \ 

k  (lO  +  S1+)  (to -  Q.)  ' 


are  included. 
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2.5  Solution  for  the  Second  Order  Field  and  Ion  Energy  Absorption 


Equation  (2d5)  for  the  second  order  fields  is 


ffUr 

tO 


(2.107) 


Through  the  use  of  equations  (2.99),  (2.100),  (2.106),  and  (2.107), 
the  relationship  between  the  second  order  electric  field  and  the 
relevant  currents  may  be  expressed  in  matrix  form  as  follows: 
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*2 

1  J  3dtj 

^31 

VC 

3  2 

II 

^32 

ii 

£  E 
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-  (i- 
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(l-od^^)  J  +  <•  V  OV  Vdx 

m  ®  c/to  K 

=  refractive  index  of  plasma  for  the  difference  frequency 
wave , 


The  various  field  components  are  then  given  by 

z 

^  z j  =  Zl  (-1)  / def  Lv] 

K=  1 


.  .  . (2.110) 


where  j  =1,  2,  3  and  denotes  components  in  the  x,  y,  and  z  directions 
respectively.  The  terms  are  the  cofactors  for  the  matrix 

[Y]  given  by  equation  (2,108), 

cle-V  tvl  =  2  (-0  VK1 

K  - 1 

The  various  cofactors  are: 
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From  equation  (2.99),  the  induced  second  order  ion  current  is 

J3  =•  tf-z  ■  ^3  ...(2.112) 

4rr!T2+- 

The  energy  absorbed  by  the  ions  from  the  second  order  electric  field, 

E^i  is  given  by 

i  +* 

W3  =  Real  (J  0  E  ) 

+&  + 

where  J_  is  the  complex  conjugate  of  .J  .  By  substituting  the  results 
from  equations  (2.100)  and  (2.112)  into  this  expression,  it  may  be 
shown  that 


=  co 
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2  2. 
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where 


H 


Similarly,  the  electron  energy  absorption  from  the  mixed  wave  may 


be  shown  to  be 


W3  =  Real  (J^  .  E3) 
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Equation  (20113)  gives  the  total  ion  energy  absorption,  that  is, 
contributions  from  both  collisionless  and  collisional  damping.  By  the 
use  of  equation  (2,100),  the  contributions  from  cyclotron  damping  to 


may  be  shown  to  be 
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where 


( ot-i)  -  ( 17+) r^T 
Kt  v  2KT+ 

(A  —  (  1  -t  A  )  S"?+ 


CHAPTER  3  CALCULATIONS  FOR  THE  PROPOSED  EXPERIMENT 


In  this  chapter,  experimental  considerations  such  as  the  selection 
of  the  D-C  magnetic  field  (Bo) ,  the  choice  for  the  incident  frequencies 
(jO^  and  oo^,  and  the  resulting  sensitivities  to  frequency,  density, 
and  angular  disturbances  will  be  discussed.  Magnetized  plasmas  with 


3.1  Selection  of  the  D-C  Magnetic  Field 

The  incident  waves  in  the  plasma  will  be  only  slightly  Landau 
damped  by  the  electrons  if  the  phase  velocities  of  the  incident  waves 
in  the  direction  of  the  static  magnetic  field  are  much  greater  than 
the  electron  thermal  velocity.  The  role  of  the  collisionless  damping 
effects  on  the  incident  waves  may  be  examined  more  closely  by  writing 
equation  (2.78)  for  the  incident  waves. 


(3.1) 


where 


(OJ;  a*3  [ 


i  =  1,2  and  denotes  the  first  and  second  incident  waves 


respectively. 


In  the  evaluation  of  the  dispersion  relations  for  the  incident 
waves,  mobility  tensors  similar  to  that  given  by  equation  (2.99)  for 
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the  difference  frequency  wave,  may  be  definedo  The  various  components 


of  this  tensor  will  be  expressed  in  terms  of  FqS  F  ,  (F1 ,  F2„ 

may  be  expressed  in  terms  of  F  by  the  use  of  equation  ( 2 0  69)  and  (2,70)), 
The  first  term  on  the  RHS  of  equation  (3d)  gives  the  collisionless 
damping  effects,  while  the  first  component  in  the  second  term  gives  the 
normal  cold  plasma  contribution,,  Only  the  latter  term  was  used  in 
evaluating  the  mobility  tensor  for  the  incident  waves  in  Section  (2,20.1) 
(see  equation  (2,12)),  It  may  be  shown  from  equation  (3d),  that  the 
effects  of  collisionless  damping  on  the  incident  waves  will  be  small 
if  |a^  J  >:>  1°  The  effects  of  electron  Landau  damping  are  represented 
by  n  -  0,  where  °v  ,  _  .  is  just  the  ratio  of  the  phase  velocity  of  the 
incident  wave  in  the  plasma  to  the  electron  thermal  velocity „  The 
case  where  n  ^  0  represents  the  contribution  of  the  n’th  harmonic  of 
cyclotron  damping. 


If  is  of  the  order  of  3,  then  only  about  one  particle  In 

2400  will  have  a  velocity  greater  than  or  equal  to  the  phase  velocity  of 


the  incident  wave 


27 


By  choosing 


+ 


a 


i(n) 


>  5,  the  effects  of  collislon- 


less  damping  on  the  incident  waves  can  be  regarded  as  negligible.  There- 
fore  the  static  magnetic  field  must  be  selected  such  that  the  dispersion 
relations  for  the  incident  waves  In  the  plasma  satisfy  the  condition 


that 


a 


i(n)l 


5,  This  would  then  insure  that  the  Landau  and  cyclo¬ 


tron  damping  of  the  incident  waves  is  negligible. 

In  order  that  the  static  magnetic  field  (B  )  may  be  selected 


such  that 


a 


>  5 9  an  estimate  for  the  maximum  electron  temperature 


i(n)l  ^ 

that  can  be  realized  in  an  experiment  is  required.  This  will  now  be 


calculated. 
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If  ion  effect  are  neglected,  the  ratio  of  the  collisionless 
to  the  collisional  damping  of  the  incident  waves  may  be  shown  to  be 
equal  to 


(3.2) 


For  the  case  in  which  [a 


_  lc  .  3/e 

^  3,  9  X  10  (  0);  +  onQ.  )  T 


collisionless  damping 

collisional  damping 


The  above  ratio  does  not  exceed  one  for  plasmas  with  densitites  greater 


the  frequencies  of  the  incident  waves  are  below  the  infrared  region . 


from  the  incident  fields  essentially  through  collisional  damping.  The 
collisional  heating  of  electrons  by  the  second  order  fields  will  be 
neglected  since  these  fields  are  much  smaller  than  the  electric  fields 
associated  with  the  incident  waves  in  the  cases  considered.  The  maximum 
electron  temperature  that  can  be  realized  in  an  experiment  is  bounded  by 
one  of  two  factors.  These  are: 

(i)  Since  the  dependence  of  the  electron-ion  collision  frequency 

-3/2 

upon  temperature  is  proportional  to  T  ,  the  rate  at  which  elec¬ 
trons  are  heated  by  the  incident  fields  will  decrease  with 
increasing  temperature.  The  maximum  electron  temperature  that  can 
be  realized  in  an  experiment  will  be  that  for  which  the  energy 
radiated  by  the  electrons  through  Bremsstrahlung  and  cyclotron 
radiation  equals  the  energy  absorbed  by  the  electrons  through  the 
collisional  damping  of  the  incident  waves. 
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(ii)  The  maximum  electron  temperature  that  can  be  realized  in 
an  experiment  will  also  be  bounded  by  the  finite,  time  that  a 
plasma  may  be  confined .  That  is,  if  the  plasma  is  confined  for 
only  one  mi  Hi- second,  the  electron  temperature  may  not  reach  the 
limit  given  by  condition  (i)  above,  where  the  rate  of  energy 
absorption  by  the  electrons  through  the  collisional  damping  of 
the  incident  waves  equals  the  energy  radiated  through  bremsstrah- 
lung  and  cyclotron  radiation „ 

By  performing  a  cold  plasma  analysis  on  the  incident  waves,  the 
power  absorbed  by  the  electrons  from  the  incident  waves  may  be  shown 
to  be 

,  e  r  s  P  ,  1 2 

W1  -  A  1  fu)j-S2.) _ Fik+,’£iu 

A  (T  l  'Si.') S  (lA±-  (b i 

.  » . (3,3) 

^  &■+) _  !Flx~  <  E±y|  +  — ^2  l^izj  ]  /Sac. 

2  8 

v  „  =  electon-ion  collision  frequency  ^  31.7no  (see  Clavier  ) 

t3/2 

Since  only  a  RH  elliptically  polarized  wave  is  transmitted  into  the 
plasma,  the  expression  for  given  by  equation  (3.3)  may  be  approximated 
as  follows,  if  terms  of  the  order  (ft  /u^)  are  neglected: 

g 

W®  —  tJpf  /  |E\J  /cm2  /sec,  „  o  o  (3<,4) 

- 

rr 


If  for  the  moment,  the  effects  of  bremsstrahlung  and  cyclotron 


-84- 


radiation  are  neglected,  the  resulting  ordinary  differential  equation 
governing  the  electron  temperature  during  the  heating  cycle  is  given  by 


/  _3_  n  \<t)  ^  40,3  u)*_  no  Ifc 

d+  ^  2  /  Cr')3/e 


...(3.5) 


where  the  energy  gained  from  each  of  the  two  incident  waves  is  taken 
to  be  approximately  equal.  In  writing  the  above  equation,  it  is 
assumed  that  the  heating  that  occurs  in  the  time,  2tt /  >  is  small. 


Therefore,  the  time  average  of  the  electron  heating  by  the  incident 
waves  appears  on  the  RHS  of  equation  (3.5).  If  the  initial  plasma 
temperature  is  approximated  as  zero, 


_ _ x57a  _  i7  ,  g  ,s 

(T$)  ~  4.&7  x  lo  u p_  x  he«.+  > 

(td±-  S2-) 


V\^ 


period  (  Sec.') 


.  .  ,  (3.6) 


where  =  electron  temperature  at  the  end  of  the  heating  period. 

2  2 

An  estimate  for  the  value  of  w  ^_/(w^  -  fl_)  is  required  if  equation 
(3.6)  is  to  be  used  for  finding  the  maximum  electron  temperature. 


i(o) 


This  information  may  be  obtained  from  the  conditions  that  a 
must  be  greater  than  or  equal  to  five  if  the  electron  Landau  and 
cyclotron  damping  of  the  incident  waves  in  the  plasma  is  to  be 
negligible.  The  ratio  of  the  two  parameters  is  given  by 


and 
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where  w  =  C  9,  =  C  , 

L  L  H  1  \ 


.Q 

+  — 


If  the  condition  that  both  a,,  A  and  a  ,  .  J  be  greater  than  five  is  to 

I  l(o)"  I  i(-l)l 

be  satisfied,  the  static  magnetic  field  (B  )  must  be  chosen  such  that 

o 

the  following  conditions  are  satisfied,, 


cU 


i(-T> 


i  lo) 


>  S’ 


>  s 


IOV 


■oY' 


i"  u)  j  ^  ^7  _ 
2  1 


^  £  Cl* 


. (3.8) 


For  instance,  for  <  w^  <  fi_/2,  if  a/i  (Q)  —  5,  then  by  equation  (3.7) 

a.,  N  /  aw  -v  <  1,  so  that  awix  is  also  greater  than  or  equal  to  5, 
l(o)  1(-1)  1(1) 

The  use  of  the  dispersion  relation  for  the  first  incident  wave 
in  equation  (3.8)  yields 


(i)  fl+  £  iOd  6.  fl-  /s. 


(0 


E= 


£<  (C  +-  (i  -  l,19xlo 


(cdi-Q.y 


{c<  -  42.  £>) 


S 


T‘ 


(ii)  CU  £  to,  ^  Q_ 

2 


...(3.9) 


a  \  /  x£  a 

(bp_  £  Cc±  +  1  /42. 8)  I  1  -  (Cj  -4g.&)  1.13  X  lo 

(cbi-Q.)2  CCi  -42.8*)  '  dq2  T  max. 


The  bound  on  maximum  electron  temperature,  resulting  from  the  finite 

time  a  plasma  may  be  confined,  is  obtained  by  using  the  results  of 
» 

equation  (3.9)  in  equation  (3.6). 

_  g 

If  terms  of  the  order  (T  /1.2  x  10  )  are  neglected,  the 

max 


following  estimates  for  the  maximum  electron  temperature  may  be 


made . 
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(i)  <  S2./2. 


T  z 

Max 


3,  64  x  10 

2“S 


Sb 


C  C±  +  o.  02 z) 
(  42.8  - 


£:. 


line. 


2 

X  heal  tn  3 


pevi  od 


(See) 


e/7 


°  Kelyia. 


o  a  o  (3  O  10) 


(ii)  ft-/e  ^  u>1  <  S2- 


T 


Wftx 


3,  x  |pgfe  (C1+  o.  023)^2. 
2f 


|e;mc|  x  hea+i^  period  (s<rc) 


2 H 


c 4 


1  k’eluiVi 


where  E  -  maximum  value  of  the  incident  electric  field  in  stat- 
volts/em.  (Note  that  the  Incident  fields  are  assumed 
to  be  RH  circularly  polarized) 0 

The  bound  on  the  maximum  electron  temperature  resulting  from  a 

balance  between  the  energy  absorbed  by  the  electrons  from  the  incident 

waves  and  the  energy  radiated  through  bremsstrahlung  and  cyclotron 

radiation  at  the  maximum  electron  temperature  will  now  be  developed. 

The  relations  for  the  energy  radiated  through  bremsstrahlung  and  cyclo- 

3 

tron  radiation  are  taken  from  Rose  and  Clark  „ 

-27  £  -  1/2  3 

W,  =  bremsstrahlung  radiation  -  1.41  x  10  Z  H.  n  (T  )  '  "  ergs/cm  /sec 

b  ~ 

. o  . (3.11) 

W  =  cyclotron  radiation  -  5.35  x  10  n_  T  (1  +  T  +. . . ) 

c  0  10 

2.4  x  10 

ergs /cm^/ sec. 

where  B  =  magnetic  field  in  gauss 

3 

n.  =  ion  density /cm 
Hr 

3 

n  -  electron  density/cm" 

T~  =  electron  temperature  in  degrees  Kelvin. 


■ 


-87- 


For  the  large  values  of  static  magnetic  field  being  considered, 
electron  cyclotron  radiation  will  dominate  over  bremsstrahlung 
radiation.  A  balance  is  obtained  by  equating  to  W^.  With  the 
help  of  equation  (3.9),  this  gives  the  following  limits  on  the 
maximum  electron  temperature: 


(i)  C2+  <  u>±  £  Q~ /£ 


I'/lAX 


1.59  X  10 


44 


C 


i 


(C)‘ 


'K 


£lvm 


(ii) 


Q-  £  ^  Sl_ 

a  1 


. . . (3.12) 


T"  a  £ 

I MAX 


44 


1.59  x  ID  (  C±  -  42 


jO  (e 


v  s^l 
i  nc-J 


C± 


M0 


2/9 


0  Kel  vm 


In  deriving  the  expressions  for  T  in  equations  (3.10)  and 

(3.12),  the  reflection  of  the  incident  waves  at  the  plasma-vacuum 

interface  has  been  neglected.  This  would  imply,  that  in  an  experiment, 

the  actual  value  of  T  will  be  slightly  less  than  that  given  by 

max 

these  equations.  This  will  result  in  the  value  for  the  D-C  magnetic 
field  (Bq)  being  slightly  larger  than  is  required  for  satisfying 
the  conditions  given  by  equation  (3.8). 

Equations  (3.9),  (3.10)  and  (3.12)  may  now  be  used  to  give  the 
required  value  for  Bq,  so  that  the  collisionless  damping  of  the 
incident  waves  in  the  plasma  is  negligible. 


' 
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(i)  S2+  <  <  Q.-/Z 


2 

Bn  > 


_  e 

as  -*  io  n0 


(  C*  +  :Ar^(c±~  z,e'e)  (l  ~  ^^4  ) 

5*  T i 


m  x 


(3.13) 


(ii)  Q- / £  -  cOi  <  S^- 


2 

Bo 


i.  8a  x  /o  up 


( +  ^q')  (^1-  42.  a)  r  1  -  (c±-  *iz.a)  2.Q7  X  io?  I 

^  ci2 


The  value  for  T  that  is  to  be  substituted  into  the  expressions 

max  r 

2 

for  B  in  equation  (3.13)  is  the  minimum  of  the  values  obtained  for 
o 


max 


in  equations  (3.10)  and  (3.12).  It  should  be  noted  that  the 


equation  for  B  when  Q  /2  <  u  <  ft  ,  places  an  upper  bound  on  an  , 
o  —  —  i  —  1 

2  9 

namely  that  (c^  -  42.8)  2.97  x  .10  >  1.  That  is,  for  sufficiently 


25  T 


max 


close  to  it  will  not  be  possible  to  satisfy  the  condition  that 
la.l(  1)1—  ^  ky  increasing  B  ,  since  the  phase  velocities  for  the  incident 
waves  in  the  plasma  are  less  than  or  equal  to  the  speed  of  light. 


3.2  Techniques  for  Optimizing  Ion  energy  Absorption 


Three  techniques  are  available  for  optimizing  the  energy 
absorbed  by  the  Ions  from  the  second  order  fields.  These  are 

(i)  By  allowing  the  difference  frequency,  w,  to  approach  the 
ion  cyclotron  frequency,  a  resonance  in  the  second  order  ion 
current  is  obtained.  The  peak  in  this  resonance  is  bounded  by 


. 


■ 
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either  collisional  or  cyclotron  damping  effects,  depending  upon 
which  process  is  dominant .  This  technique  places  restrictions  on 
the  allowable  frequency  and  D-C  magnetic  field  fluctuations.  The 
.  effect  of  density  fluctuations  is  small. 

(ii)  Another  technique  for  optimizing  the  ion  energy  absorption 
from  the  second  order  fields  is  to  obtain  a  resonance  in  the 
magnitude  of  the  second  order  fields ,  by  allowing  the  mixed  wave 
to  approach  a  natural  mode  in  the  plasma.  If  to  is  sufficiently 
removed  from  the  ion  cyclotron  frequency,  the  resonance  will  be 
limited  by  collisional  effects.  This  method  however  imposes 
stringent  restrictions  on  the  allowable  density  fluctuations  and 
the  allowable  perturbations  in  the  angle  of  incidence  for  the 
sources  (see  Figure  (2.2)),  in  addition  to  restrictions  on 
frequency  and  static  magnetic  field  disturbances. 

(iii)  A  third  possible  scheme  is  the  combination  of  the  above  two 
techniques.  However,  because  a  field  resonance  effect  is  utilized, 
this  method  of  ion  heating  will  be  very  sensitive  to  density, 
frequency,  static  magnetic  field,  and  angular  disturbances. 

Examples  of  the  first  and  third  schemes  will  be  given  in  section 

(3.4).  Examples  of  the  second  scheme  have  been  given  by  James  and 
17  L  H 

Thompson  ,  and  Jayasimha'  „  These  will  be  discussed  in  Chapter  4. 


3.3  Anisotropic  Pressure  Effects 


When  the  first  and  third  schemes  in  Section  (3.2)  are  used, 


-90- 


the  energy  from  the  mixed  wave  is  selectively  coupled  to  the  ions 
in  a  direction  perpendicular  to  the  static  magnetic  field.  Rapid 
randomization  of  the  ordered  motion  perpendicular  to  the  magnetic 

field  will  then  occur  through  fine  scale  mixing  as  described  by 

1  £9 

Berger  et  al  and  Stix  ,  increasing  the  temperature  of  the  trans¬ 
verse  motion  of  the  ions.  That  is,  the  thermal  motion  of  the  ions  in 
the  z-direction  will  cause  ions  which  are  in  one  region  of  the  wave  to 
move  into  another  region  where  the  phase  of  the  wave  is  different. 

A  large  number  of  such  transitions  will  tend  to  increase  the  temper¬ 
ature  of  the  transverse  ion  motion.  Once  the  phase  mixing  has 
occured,  collisional  effects  will  tend  to  return  the  plasma  to 
Maxwellian.  However,  the  maximum  perpendicular  ion  temperature  (T^) 

that  may  be  realized  is  bounded  by  the  ion  temperature  parallel  to 

+ 

the  magnetic  field  (Tft )  by  the  onset  of  an  ion  cyclotron  overstability. 
The  limit  on  T*  is  approximately  given  by 


± 


<  (3 


(3.14) 


T. 


where 


(2  =  e-rr  Ki  K' 


ion  kinetic  pressure  parallel  to  B 


Bo 


magnetic  pressure 


o 


For  the  examples  to  be  considered  in  the  next  section,  the  ratio  of 

+  +  cannot  exceed  twenty.  This  is  however  a  pessimistic  estimate 
Tj.  /Tu 

for  the  maximum  ion  temperature  because  the  ion  temperature  parallel 
to  the  magnetic  field  will  increase  during  the  heating  cycle  through 
collisional  effects.  To  extend  the  upper  limit  on  the  ion  temperature 
as  given  by  equation  (3.14),  the  plasma  may  be  heated  by  a  sequence 
of  pulses.  The  period  between  the  pulses  could  then  be  arranged  to  be 


-eneiJ  9/13  Jo  97ii3oJ3qs93  sriJ  gnieesu/ii  ,  '*i3‘  boi.  J  >-  3  9  703798 
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of  sufficient  length  so  that  the  ions  reach  equilibrium  before  each 
heating  pulse,, 

The  power  absorbed  by  the  ions  through  the  collisionless  damping 
of  the  second  order  fields  for  the  case  where  Ta  ^  Tu  may  be  shown 
to  be 


I  g 

W,  (  Colli  sunless)  =  cOp+  \<T1 


2“  Jrr  K-z,  wi+ 


3/2 


(co-Qv)  ix  +  h+tu 

to  ~t*  i 


,  +  2 
enp-  (  oC_ 


(3.15) 


+  1 


£  2 

4  (u>-  H-t-')  /  k*  ^  |  h 

51+2  ' 


3Z 


4  £(<0- 


— ~)  il*  i-re  ^  +  1  ))  | 

S2+-  j 


Equation  (2.115)  may  be  obtained  by  setting  Tj_  -  Tu  in  equation  (3,15). 

+ 

If  Tn  is  taken  to  be  constant,  and  if  a>  is  sufficiently  close  to 

it  may  be  seen  from  equation  (3.15)  that  the  power  absorbed  by  the 

ions  through  the  cyclotron  damping  of  the  difference  frequency  wave  will 

“f” 

remain  essentially  constant  with  Increasing  Tx<,  From  the  above  equation 
it  may  also  be  seen  that  the  ions  will  absorb  energy  from  the  driven 
second  order  fields  only  if 


( co-  fl+A  Tj_  4-  £2  +  ~T\\ 


>  O 


. . . (3.16) 


This  requires  that 


e  >  -r*/ T+ 


where  tl>  -  (l4A).Q+ 


By  using  the  limit  established  by  the  condition  of  a  cyclotron 
overstability  for  T^/T^  (equation  (3.14),  the  restriction  on 
A  in  equation  (3.16)  becomes 


A  > 


o  o 


. (3.17) 
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3,4  Examples 


In  this  section,  the  theory  developed  in  Chapter  2  and  in 

Sections  (3.1)  and  (3,2)  of  this  chapter,  will  be  applied  to  plasmas 

12  3  15  3 

with  densities  in  the  range  10  /cm  to  10  /cm  ,  For  the  cases 
considered,  the  square  of  the  ratio  of  the  ion  Larmor  radius  to  the 
perpendicular  wavelength  will  have  a  value  less  than  0.5  for  a 

g 

maximum  ion  temperature  of  10  degrees  Kelvin.  This  ratio  will  be 
represented  by  X+.  In  the  examples  considered,  the  incident  wave 
outside  the  plasma  will  be  taken  to  be  RH  circularly  polarized,  with 
an  electric  field  of  magnitude  of  eY  .  This  paricular  wave  is 
chosen  for  better  matching  at  the  plasma-vacuum  interface,  since  the 
transmitted  wave  is  almost  RH  circularly  polarized.  The  numerical 
calculations  were  performed  on  an  IBM  360  computer. 

In  Figure  (3.1),  the  static  magnetic  field  (B^)  and  the  value 

i  s 

for  maximized  with  respect  to  the  angle  of  incidence  9^  (see 

Figure  (2.2))  for  each  without  the  use  of  field  resonance  effects, 


have  been  plotted  as  a  function  of  C^,  where  co^  = 

It  is  evident  from  this  Figure  that  a  good  choice  for  the  frequency 
of  the  incident  waves  is  in  the  neighborhood  of  the  hybrid  frequency. 

If  the  static  magnetic  field  is  at  a  premium,  a  higher  incident 
frequency  may  be  used.  However,  unless  a  field  resonance  effect  is 
utilized,  this  will  result  in  a  reduction  in  the  rate  of  ion  energy 
absorption. 

Figure  (3,2)  shows  the  dependence  of  on  the  angle  of  incidence 
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for  the  first  source  (see  Figure  (2.2)),  with  u)^  being  equal  to 

1-1  tunes  the  hybrid  frequency,,  The  sharp  peak  in  for 

^  35  is  the  result  of  a  field  resonance  effect .  This  resonance  is 

strongly  dependent  upon  density,  angular,  frequency  and  magnetic 

field  fluctuations  as  is  shown  by  Figures  (3.3),  (3.4),  and  (3.5) , 

If  the  sensitivities  are  defined  using  the  one-quarter  power  points 

17 

as  in  James  and  Thompson  9  then  the  following  sensitivities  may  be 
defined  for  the  above  field  resonance, 

6to^/oo^  .67  x  10 

6n  /n  ^  7  x  10  ^ 
o  o 

6B  /B  v  3d  x  10“ 4  .,.(3.18) 

O  o 

6(0^)  ^  10  4  radians 

for  T  =  106  °K 


The  dependence  of  on  the  detuning  of  the  difference  frequency 
from  the  ion  cyclotron  frequency  and  the  effects  of  density  and  static 
magnetic  field  perturbations  is  given  in  Figures  (3.6)  and  (3.7)  for 
0^  =  64 o 1° .  If  the  ions  are  to  gain  energy  through  cyclotron  damping, 
the  above  choice  for  an  operating  point  represents  the  worst  case  for 
sensitivity  to  frequency  and  static  magnetic  field  fluctuations. 

The  value  for  in  this  case  is  just  large  enough  so  that  cyclotron 


“3  l 

damping  occurs  for  (co  -  fi+)/fi4_  ^  10  .  That  is  |a 


1  +  1 
a  ,  = 

(a)  -  ft+) 

m 

1  -l1 

K  \ 

+ 

Z  1 

2KT- 

for  (w  -  )/fi  r ^  10 

*r  "r 


-5 


If  w  is  detuned  any  further  from  ft  ,  the 


'v  1 


difference  frequency  wave  will  not  be  cyclotron  damped.  This  is  the 
reason  for  the  high  sensitivity  to  frequency  and  static  magnetic  field 
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fluc.tuations  at  this  operating  point, 

A  much  better  choice  for  an  operating  point  is  given  by 
s 

=  50° 2  ,  The  results  for  this  case  are  given  by  Figures  (3,8) 

and  (3,9),  In  this  case,  the  absolute  value  of  K  is  much  larger 

z 

s 

than  that  obtained  when  ©^  =  64,1''0  This  may  be  seen  by  examining 


the  values  for  la  I  (a 

‘O'  o 


0J 


/K 


m 


)  in  Figure  (3,10)  for  the  case 


2KT" 


where  the  slope  of  the  plasma-vacuum  interface  is  zero,  and  ©~  has 
the  values  50,2'  and  64, lc,  This  would  imply  that  for  the  operating 
point  given  by  ©^  =  50,2°,  the  resulting  sensitivities  to  frequency 
and  static  magnetic  field  fluctuations  will  not  be  as  critical  as 

S 

in  the  case  where  ©^  =  64,1°,  The  sensitivities  obtained  from  the 
quarter  power  points  are  as  follows  for  the  case  where  ©^  =  50,2°, 


ScOf/^i  ^  1,2  x  10 

^  n  /n  >  ,  2 
o  o 

/B  %  5,6  x  10“ 4 
o  o 

6(©^)  <  30°  radians 
for  T  =  106  °K 


(3,19) 


By  using  the  results  in  Figure  (3,10),  it  may  be  shown  that  plus  and 
minus  five  degree  fluctuations  in  the  slope  of  the  plasma  vacuum 
interface  may  be  tolerated  without  leading  to  excessive  electron  Landau 
damping  of  the  mixed  wave.  From  Figure  (3,6)  and  (3,8),  it  is  evident 
that  with  increasing  temperature,  the  rate  at  which  energy  is  absorbed 
by  the  ions  will  drop  slightly,  however  an  improvement  will  be  obtained 


. 
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in  the  frequency  bandwidth  over  which  the  mixed  wave  is  cyclotron 
damped „ 

The  field  resonance  for  to1  =  lOfi  is  investigated  in  Figures 

-L  H 

(3,12)  to  (3,14),  If  the  sensitivities  are  defined  by  the  one- 
quarter  power  points  as  in  equation  (3,11),  then  for  the  above 
resonance 


^  9.1  x  10 

6n  /n  '’u  7,4  x  10  ^ 
o  o 

6B  /B  ^  3o 9  x  10~4  . . . (3.20) 

o  o 

6(0^)  _<  1.6  x  10  ^  radians 

for  T  -  106  °K 


The  field  resonance  for  oa  -  10ft  is  considerably  "sharper"  than 

.1  H 

that  for  oo1  =  hlfl  ,  resulting  in  tighter  restrictions  for  the 
.1  H 


allowable  frequency  and  angular  disturbances  on  the  incident  waves, 

.15  3 

The  results  for  a  plasma  with  a  density  of  10'  /cm  is  given 


in  Figures  (3.15)  to  (3,21). 
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Source.  The  physical  parameters  are  as  given  in  Figure  (3.2). 
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negative  6n  and  6B  are  given  by  the  broken  curves. 
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given  by  the  broken 
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negative  6n  and  6B  are  given  by  the  broken 
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are  given  by  the  broken 
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negative  6n  and  6B  are  given  by  the  broken 


-105- 


(3.1),  with  C  =  1.1.  The  values  of  9  are  given  for  the  condition  that 
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for  the  first  source.  The  physical  parameters  are  as  given  in  Figure  (3.11). 
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CHAPTER  4  SUMMARY  AND  DISCUSSION  OF  HEATING  TECHNIQUES  USING  THE 


NONLINEAR  MIXING  OF  TWO  WAVES 


In  this  chapter,  a  summary  of  the  results  pertaining  to  ion 
heating  through  the  non-linear  mixing  of  two  waves  will  be  given. 

The  merits  of  the  techniques  that  may  be  used  to  optimize  the  energy 
absorbed  by  the  ions  from  the  second  order  fields  will  be  discussed. 


4.1  HF  Heating  of  Ions  by  Maximizing  the  Second  Order  Fields 


Two  schemes  have  been  suggested  recently  (James  and  Thompson  , 

18 

and  Jayasimha  )  where  the  collisional  heating  of  ions  is  optimized  by 
maximizing  the  second  order  fields.  A  resonance  in  the  magnitude  of 

the  second  order  fields  is  obtained  by  allowing  the  mixed  wave  to 

approach  a  natural  mode  in  the  plasma,  which  in  the  above  case  is  an 

extra-ordinary  Alfven  wave  propagating  perpendicularly  to  the  static 

magnetic  field. 


4.1.1  HF  Heating  of  Ions 


Analysis  by  James  and  Thompson"7  ) 


The  field  vector  diagram  for  the  ion  heating  scheme  suggested  by 
James  and  Thompson"7  is  given  in  Figure  (4.1). 
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FIGURE  4.1.  Field  Vector  Diagram  for  HF  Ion  Heating.  This  is  essentially 
Figure  (1)  in  James  and  Thompson  ^  however  the  exponential  dep¬ 
endence  of  the  field  quantities  is  written  as  to  agree  with  the 
notation  used  in  equation  (2.3).  The  plasma  is  assumed  to  occupy 
the  half  space  z  >0. 


The  authors  used  a  cold  plasma  model  to  describe  the  nonlinear 
generation  of  the  second  order  fields  and  the  resulting  ion  heating. 
An  investigation  into  the  sensitivity  of  the  heating  process  to 
fluctuations  in  frequency,  density,  and  d-c  magnetic  field  was 
carried  out.  These  results  will  be  compared  with  those  obtained  by 
the  use  of  Whistlers  in  Table  (4.1).  The  sensitivity  of  the  heating 
scheme  suggested  by  James  and  Thompson"^  to  fluctuations  in  the  angle 
of  incidence  for  the  first  wave  (0  in  Figure  (4.1))  will  now  be 
calculated.  The  propagation  constant  for  an  extra-ordinary  Alfven 
wave  propagating  perpendicularly  to  the  static  magnetic  field  is 
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Equation  (4.1)  and  the  requirement  that  the  mixed  wave  propagate 
perpendicularly  to  the  static  magnetic  field  may  be  used  to  give 
the  required  value  for  the  difference  frequency  as^ 
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If  9j_  is  the  required  angle  of  incidence  for  the  first  wave  so  that  K 
is  perpendicular  to  the  static  magnetic  field,  the  use  of  the  free 
space  values  for  the  propagation  vectors  and  in  Figure  (4.1) 
yields 
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If  the  error  in  setting  the  required  value  of  0  is  0  ,  then  for 
small  values  of  0  ,  the  x-component  of  the  propagation  vector  for 
the  mixed  wave  is 

K'x  =  ^t/c  Sm  (  0j_  f 


/c  Sm  0J. 
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1  Alfven 
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No  resonance  effects  will  be  observed  when 

&C  /  ©j_  01 1) 

Equations  (4*2)  and  (4o3)  may  be  used  to  show  that  the  above  condition 
is  equivalent  to 


©  £  Q.+  SI  -  /id  p= 


(4.5) 


If  the  sensitivity  of  the  heating  process  to  angular  fluctuations  is 
defined  by  the  quarter-power  points,  the  analysis  given  by  James  and 
Thompson"^  may  be  extended  to  give 
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Condition  (4o6)  is  much  more  severe  than  is  (4.5).  For  the  example 
considered  in  [17],  (4,6)  gives 

0r  icfa  radons 

C  ...(4.7) 

For  the  same  case,  if  any  resonance  effect  is  to  be  observed  at  all, 
equation  (405)  gives  the  requirement  that 

00  «•  Ob  Tadiuns  C  ^  d  eg  Y~e  e.%)  ,,  o\ 
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4 o.l .  2  HF  Heating  of  Ions  (Analysis  by  Jayasimha  ) 


In  this  analysis,  two  incident  waves  with  frequencies  approxim¬ 
ately  equal  to  the  electron  plasma  frequency  are  sent  into  the  plasma 
at  right  angles  to  the  static  magnetic  field.  The  geometry  of  the  plasma 
and  the  sources  is  given  in  Figure  (4.2) 


Y 


Figure  4,2.  Field  Vector  Diagram  for  HF  ion  Heating  used  by  Jayasimha 
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The  results  obtained  for  this  case  are  closely  related  to  those 

17 

obtained  by  James  and  Thompson  ' „  By  the  use  of  lower  frequency 

incident  waves  ,  a  gain  of  approximately  four  hundred  is  claimed  in 

the  power  absorbed  by  the  ions  from  the  mixed  wave.  An  improvement  in 

the  allowable  fluctuations  in  the  angular  position  of  the  sources  is 

also  realized  through  the  use  of  lower  frequency  waves.  However,  this 

gain  in  the  rate  of  energy  absorption  by  the  ions  may  not  be  realized 

in  practice  because  of  the  tighter  restrictions  on  frequency,  density, 

and  d-c  magnetic  field  fluctuations „  The  sensitivity  of  the  heating 

process  to  frequency,  density  and  d-c  magnetic  field  fluctuations  are 

12  3 

summarized  below  for  a  plasma  with  a  density  of  10  /  cm  and  a  temp¬ 

erature  of  10^  degrees  Kelvin „ 


^  rA-  Kl 

—  -3.4*10 

(a)  2. 

M  4*  Cc)p_ 

‘s’t'lo 

2nJ  Ae.\ 

_  fo 

10 

Ho 

S’ Bo 

Kl  As' 

-  <o 

^  I  o 

Bo 

S2+ 

where  N  =  w/fl  .  The  above  results  are  an  order  of  magnitude  more 

sensitive  to  frequency,  density,  and  d-c  magnetic  field  fluctuations 

17 

than  those  obtained  by  James  and  Thompson 

The  finite  Larmor  radius  effects  would  appear  to  play  a  role  more 
important  than  is  suggested  because  the  dominant  electron  driving  terms 
cancel,  leaving  terms  of  order  (m_/m  ) ,  and  this  is  of  the  same  order 
as  the  dominant  ion  driving  terms. 


402  Comparison  Between  HF  and  Whistler  Heating 


The  use  of  high  frequency  waves  (go  ,0)2  >  )  as  suggested 

17 

by  James  and  Thompson  is  compared  to  the  use  of  Whistler  waves 

<  a)  <  9,_)  as  suggested  in  this  thesis,  in  Table  (4.1).  A 

comparison  is  made  between  the  results  obtained  for  heating  plasmas 

12  3  6 

with  a  density  of  10  /cm'  and  a  temperature  of  10  “Kelvin.  Two 

possible  modes  of  operation  for  Whistler  heating  are  presented „  One 

is  the  use  of  the  resonance  obtained  in  the  second  order  ion 

current  by  allowing  the  difference  frequency  to  approach  the  ion 

s 

cyclotron  frequency,.  The  physical  parameters  (n  9'^,.,) 

are  chosen  to  be  those  given  in  Figures  (3.6)  and  (3,8).  The  second 
method  uses  an  ion  current  resonance  coupled  with  the  field  resonance 
that  is  obtained  by  having  the  difference  frequency  wave  approach  a 
natural  mode  in  the  plasma „  The  physical  parameters  are  chosen  to 

1 8 

be  those  given  in  Figure  (3ol2)„  The  claims  suggested  by  Jayasimha 

17 

have  already  been  compared  to  the  results  given  by  James  and  Thompson 
in  Section  (4. 1.2)  and  will  not  be  repeated  here. 

The  results  obtained  for  the  cases  where  a  second  order  current 
and  field  resonance  are  used  to  optimize  the  power  absorbed  by  the  ions 
(eg,  see  Figure  (3,11))  will  be  valid  for  incident  electric  field 
amplitudes  less  than  3  KV/cm„  The  field  resonance  effect  will  cause 
the  second  order  fields  to  approach  the  same  order  of  magnitude  as 
the  incident  fields  for  larger  values  of  j  |  ,  giving  rise  to 
secular  terms.  This  may  be  resolved  by  considering  the  mixed  fields 
to  be  of  the  same  order  as  the  incident  fields  in  the  perturbation 
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TABLE  4.1  COMPARISON  BETWEEN  HF  AND  WHISTLER  HEATING 
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theory  „  The  perturbation  theory  used  in  this  thesis  is  valid 

for  field  intensities  to  70  KV/cm  in  the  cases  where  the  energy 
absorbed  by  the  ions  is  optimized  only  through  an  ion  current 
resonance. 


4.3  Conclusions 


From  Table  (4.1)  it  is  evident  that  the  sensitivity  of  a  heating 
process  employing  the  nonlinear  mixing  of  whistlers  to  frequency, 
density,  d-c  magnetic  field  and  angular  fluctuations  is  considerably 
less  than  that  in  which  HF  waves  are  used.  If  only  an  ion  current 
resonance  effect  is  utilized,  the  significant  gain  is  in  the  relax¬ 
ation  of  the  stringent  restrictions  on  density  and  angular  fluctu¬ 
ations.  Through  the  use  of  whistlers,  the  sensitivity  to  density 
fluctuations  is  reduced  by  four  orders  of  magnitude  while  the 
sensitivity  to  angular  fluctuations  is  reduced  by  about  seven  orders 
of  magnitude.  An  order  of  magnitude  is  gained  in  the  rate  of  ion 
energy  absorption.  The  sensitivities  to  frequency  and  d-c magnetic 
field  fluctuations  are  reduced  by  two  and  one  order  of  magnitude 
respectively.  Since  the  use  of  an  ion  current  resonance  effect 
only  requires  that  the  difference  frequency  be  sufficiently  close  to 
the  ion  cyclotron  frequency,  the  only  critical  adjustments  in  an 
experiment  will  be  that  of  the  difference  frequency  and  the  d-c 
magnetic  field. 
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If  a  field  resonance  effect  is  utilized,  in  addition  to  an  ion 
current  resonance  the  mixed  wave  must  approach  a  natural  mode  in  the 
plasma,  which  in  the  cases  considered,  is  an  extraordinary  Alfven  wave. 
This  results  in  much  tighter  restrictions  on  frequency,  density,  d-c 
magnetic  field,  and  angular  fluctuations  than  is  required  in  the  case 
where  only  an  ion  current  resonance  effect  is  utilized. 

Two  advantages  of  the  use  of  HF  waves,  rather  than  whistlers  in 
the  mixing  process  are: 

(i)  No  power  is  reflected  at  the  plasma-vacuum  interface  when 
HF  waves  are  used. 

(ii)  The  transmission  of  HF  waves  into  a  plasma  will  be  indepen¬ 
dent  of  the  plasma-vacuum  interface  geometry. 

17 

In  the  analysis  performed  by  James  and  Thompson  ,  it  was  found 

that  E„  =  ifi./to  E„  .  For  the  case  considered  in  [  17],  fi,/oo  'v-  1/8. 

3y  +  3x  + 

That  is,  the  second  order  electric  field  was  right  hand  elliptically 
polarized.  If  oo  were  allowed  to  approach  this  would  suggest  that 
the  second  order  electric  field  would  become  almost  RH  circularly 
polarized.  This  is  just  what  is  found  in  the  case  of  whistler  heating 
near  field  resonance  with  co  near  (See  Figure  (3.2).  For  both  the 
ion  cyclotron  and  collisional  damping  of  the  difference  frequency  wave, 
the  ions  resonate  only  with  the  electric  field  component  rotating  in 
the  direction  of  ion  gyration  in  the  static  magnetic  field.  This 
would  then  suggest  that  a  possible  modification  of  the  heating 
scheme  presented  in  this  thesis  would  be  to  have  the  ions  absorb 
energy  from  the  difference  frequency  wave  near  the  second  harmonic  of 
the  ion  cyclotron  frequency  so  as  to  increase  the  amount  of  LH 


■ 
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polarization  0 


This  aspect  would  certainly  warrant  further  investigation. 
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APPENDIX  A.  METHOD  OF  ORBIT  INTEGRATIONS  APPLIED  TO  SECOND  ORDER 


The  Vlasov  equation  is 


IS.  +  v  •  IS  +  (£  +  vx  s']  .  7>f  -  o 


.  .  .  (A. 1) 


Equation  (A.l)  will  be  solved  assuming  the  following  ordering  of 
variables ; 
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The  perturbed  distribution  functions  will  be  solved  for  in  a  Lagrangian 
system  of  coordinated.  The  change  of  a  distribution  function  expressed 
along  some  trajectory  defined  by  _r(t)  is 


d-P  —  1)P  +  .  dr  +  2>S  ,  d  y  (OKev"£  y  =  dr 

dt  ht  "}y  dt  ~}\/  d  t  dt 

In  this  particular  case,  the  zero'th  order  trajectory  of  a  particle 
in  a  magnetic  field  will  be  chosen,  that  is 


dr  -  V  dy_  -  zj  eCj  v  x  Bc  eA  =  sign  of  species 

dt  dt  ki  s  c 
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Substituting  (A.  2)  into  (A.l) 
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The  equations  for  the  various  perturbations  are  now  obtained  by 
using  harmonic  balance,,  For  instance,  by  equating  terms  with 
exponential  dependence  i(K0r  -  wt) , 


^£3  4  V  •  ££3  t  Z^  e€±  (  y  x  B0)  ,  c>-G  -  -  Zj  e  £]  ( t-  4  \/x  B3)  .  £fo 

~2rL  yy>  -  c  c  ~2>v 

-  £j<s£j  (£±+  M*  Bj)  .  3fe  _  Zjec^  (Eg  +  v  x  Bp)  .  Afj 
ma  C  ^  ~ ^  c  /  is 

... (A. 3) 

The  left  hand  side  of  equation  (A. 3)  is  just  the  change  of  the  distri¬ 
bution  function  f ^  along  the  zero’ th  order  trajectory  of  a  particle 
in  a  magnetic  field.  The  solution  of  equation  (A. 3)  and  its  velocity 
moments  will  give  the  second  order  current.  This  may  then  be  used 
with  Maxwells  equations  to  give  the  second  order  fields  and  the 
resulting  dissipation  can  then  be  found. 
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APPENDIX  B  METHOD  OF  SOLUTION  FOR  f 
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The  asymptotic  form  of  the  assumed  solutions  for  the  incident 
waves  is  such  that  as  t  ->  -  °°,  the  value  of  the  electric  and  magnetic 
field  components  tend  to  zero.  From  equation  (2.3),  this  implies  that 
Im(rn^)  >  0.  The  imaginary  component  of  will  be  denoted  as  un  . 

A  similar  argument  applies  to  the  second  incident  wave,  so  that  the 
imaginary  component  of  the  mixed  frequency  may  be  written  as 


uk  =  Im  (uj  )  +  ImCu^) 
From  equation  (2.10) 
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where  all  constants,  terms  linearly  dependent  upon  t1’  and  periodic 
terms  as  exp(-i  oo^t’*)  have  been  incorporated  into  the  term  denoted 
by  (  ).  The  results  given  in  equations  (2.4)  and  (2,5)  may  be 

used  to  express  the  variables  at  time  t'*  in  terms  of  their  values  at 
time  t’.  If  t'  is  defined  as  follows: 

x'  =  t'  -  t" 


equation  (B.2)  may  be  rewritten  as 
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If  the  above  expression  is  substituted  into  equation  (2.49),  then 


32 


<  i  rv\ 

T  — ->  aO 
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where  constant  and  periodic  terms  have  been  incorporated  into  the  term 
denoted  by  [  ].  The  condition  that  t' 1  <_  t '  ,  as  given  by  the 

integration  limits  in  equation  (B.2)  requires  that  >_  T  in  equation 
(B.4).  The  ambiguity  in  the  definition  of  the  upper  limit  in  the 
second  integral  in  equation  (B.4)  as  T  00  will  be  resolved  by  writing 
this  equation  as  a  sum  of  two  terms,  the  second  of  which  will  be 
shown  to  be  zero. 
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For  e  small,  but  non-zero,  and  with  the  condition  that  T^  >_  T,  the 
second  integration  in  the  first  term  in  equation  (B.5)  may  be  written 


as 


oo 
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'  o 


dr'  ( 
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As  T  ->  00 ,  the  upper  limit  tends  to  e  00  in  the  worst  case,  which  may 
be  replaced  by  00  if  eis  small  but  non-zero.  Consequently  the  first 
integration  in  the  first  term  on  the  RHS  in  equation  (B.5)  defines  an 
integration  on  the  interval  (  -°°,t].  The  second  integration  in  the 
first  term  may  be  taken  over  the  closed  interval  [  0,  00 ]  .  The  first 


. 
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term  on  the  RHS  of  equation  (B.5)  with  e  ->  0  is  that  given  by  equation 
(2.50) . 

The  absolute  value  of  the  second  term  in  (B.5)  will  be  shown  to 
be  smaller  than  any  positive  number  and  therefore  zero  as  e  -*  0.  Let 

_T7  /  t'+T, 

A  -T1£  -  I  i 

T— ■? 

T4  o*  -T 

Since  >  0,  the  absolute  value  of  the  second  integral  is  bounded. 

This  value  will  be  denoted  by  M^.  The  maximum  value  of  the  term, 

* 

E.  E.  [  ]  will  be  taken  to  be  less  than  some  finite  value  M_. 

lx  2x  2 

Then 

i'c-T)  _ cJ ; T  CO;T  \ 

~  \  dt^i^ae  =  e  (e  -1)  ...(B.7) 

J0  oJ;  v  J 

By  choosing  e  sufficiently  small,  the  RHS  of  equation  (B.7)  may  be 
made  less  than  any  positive  number.  Therefore  Af^  =  0* 
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